Zero-frequency noise in adiabatically driven, interacting quantum systems 
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We investigate current-current correlations of adiabatic charge pumping through interacting quan- 
tum dots weakly coupled to reservoirs. To calculate the zero-frequency noise for a time-dependently 
driven system, possibly in the presence of an additional stationary bias, we perform within a real- 
time diagrammatic approach a perturbative expansion in the tunnel coupling to the reservoirs. We 
apply this formalism to study the adiabatic correction to the zero-frequency noise, i.e., the pumping 
noise, in the case of a single-level quantum dot charge pump. If no stationary bias is applied, the 
adiabatic correction shows Coulomb-interaction-induced deviations from the fluctuation-dissipation 
theorem. Furthermore, we show that the adiabatic correction to the Fano factor carries informa- 
tion about the coupling asymmetry and depends of the choice of the pumping parameters. When 
including a time-dependent finite bias, we find that there can be pumping noise even if there is zero 
adiabatically pumped charge. 

PACS numbers: 72.25.-b, 73.23.Hk 



I. INTRODUCTION 

Transport through mesoscopic systems, such as quan- 
tum dots, can be achieved through a periodic time- 
dependent modulation of the system. This is referred 
to as pumping. In the regime of adiabatic pumping, 
the modulation is slow, such that the system can al- 
most immediately follow the driving. A directed trans- 
port through quantum dots in the absence of a stationary 
bias (or, alternatively, also an additional pumping con- 
tribution on top of or opposite to transport due to a 
stationary bias) is then possible only if there are at least 
two independent pumping parameters at work.ElEl 

Adiabatic charge pumping offers the possiblity of a 
controlled emission of one charge per cycle, useful as a 
quantum standard for the current or as a coherent parti- 
cle source for quantum operations ."^^^ In addition to the 
current signal, knowledge about the noise created dur- 
ing a pumping cycle is vitally important, for instance in 
metrology to investigate the quantization limitations of 
the pump ,1^^11121 and also as a tool to detect signatures for 
emission of coherent electron packets. '^^'l^ 

Recently there has also been strong interest in pump- 
ing due to quantum-interference effects, 'S^itSfiJ where the 
pumped charge is not quantized. In this context it has 
been shown that the pumped charge may be used as 
a spectroscopic tool that reveals features beyond time- 
independent measurements.'^23^ 

Another powerful tool for an extended transport spec- 
troscopy is the transport noise. In this manuscript, we 
aim to identify additional spectroscopy features in the 
pumping noise. In the time-independent case with an ex- 
ternal bias, there are two system-intrinsic noise contribu- 
tions, namely thermal (Johnson-Nyquist) and shot noise, 
where the latter can carry characteristic information 
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FIG. 1: Energy diagram of a single-level quantum dot at- 
tached to two leads via tunneling barriers. A symmetric volt- 
age bias (with respect to the equilibrium chemical potential, 
see horizontal line) can be applied across the system and the 
energy level, the bias and the coupling strength may be chosen 
time dependent. Time-dependent charging and uncharging of 
the system goes along with a time-dependent difference be- 
tween the currents flowing into the leads. 



about the nature of the charge carriers and their statis- 
tics. Various theoretical and experimental works dealt 
with the study of noise in stationary transport .'^SEH Yoy 
a time-dependently driven system an additional pumping 
noise emerges.!^ For different types of time-dependent 
driving, noise and the full counting statistics has been in- 
vestigated theoretically,I^^M ^ and has been proposed as a 
measure of quantum effects. I^SlllHIll Recently, also exper- 
iments measuring no ise in ti me-dependently driven sys- 
tems were performed! ^^ * ^^ * ^^* ^ 

When considering quantum pumps realized in small 
quantum dots. Coulomb interaction between electrons 
on the dot becomes important, due to the small ca- 
pacitance of strongly confined systems. The impact of 
Coulomb interaction on the transport behavior of adi- 
abatic pumps has recently been studied in various se- 
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tups and regimes In particular, it has been shown 
that adiabatic quantum pumps can reveal characteristic 
fingerprints of the electron-electron interaction!232ll Also 
the noise in stationary systems is sensitive to Coulomb 
interaction effects.^ The study of the noise in time- 
dependent systems with strong Coulomb interaction has 
not been addressed so far. In this paper, we aim to close 
this gap, and include Coulomb interaction with arbitrary 
strength. Based on new noise features, we extend the 
possibilities of spectroscopy and find signals that give 
evidence of the underlying pumping mechanism. 

In this paper we derive a formalism to compute the 
zero-frequency noise for adiabatically driven quantum 
systems, based on a real-time diagrammatic approach 
that enables us to include an arbitrary Coulomb interac- 
tion strength and non-equil ibrium due to a bias and due 
to time-dependent driving .E^*^ This general real-time 
formalism has previously been the basis for the develop- 
ment of a method to calculate the zero-frequency noise 
in a strongly interacting quantum dot system driven by 
a stationary biaJ^ and for the development of a method 
to calculate charge and spin transport in qu^Jntum dot 
systems with a slow time-dependent driving .'23221 

We implement the formal results of the method de- 
veloped here to study the pumping noise for the sim- 
plest generic model, a single-level quantum dot coupled 
to two leads, taking into account Coulomb interaction, 
for a wide range of possible time-dependent parameters. 
The quantum dot system is shown in Fig. [T] We are able 
to present analytic expressions for all driving regimes, al- 
lowing for a detailed analysis of the origin of the various 
noise contributions. 

We start our discussion by investigating the pump- 
ing noise when the non-equilibrium enters purely via 
the time-dependent driving. We rec over a time-averaged 
fluctuation-dissipation relatiorl^^I^ and study deviations 
from it in the presence of both interaction and time- 
dependent driving. From the correction to the Fano fac- 
tor due to pumping, we can (independent of the choice 
of pumping parameters) extract signatures of system- 
specific parameters, such as coupling asymmetry, as well 
as the level position with respect to the electron-hole 
symmetric point. 

We, furthermore, consider the case of having an addi- 
tional, possibly time-dependent, external bias. For a cer- 
tain pumping prescription we find in the low-bias regime 
that the pumping noise is finite even if the time-averaged 
pumping current, namely its dc component, vanishes. We 
show how the time-averaged (i.e., zero-frequency) pump- 
ing noise feature provides information about the time- 
resolved pumping current; namely that a finite pumping 
noise is present if there is a finite time-resolved pumping 
current, and zero pumping noise if the pumping current 
vanishes for each instant of time. 

This manuscript is structured as follows. We start by 
introducing the model and outline the derivation of the 
noise formulae in Sec. |ll] The main results, where we ap- 
ply our method to a driven single- level quantum dot, are 



divided into a discussion of pure pumping (no bias) in 
Sec. mil and the calculations for a finite bias in Sec. IIVI 
The conclusion in Sec. |V]is followed by an Appendix de- 
tailing the derivation of the pumping noise. 

A. Current and noise in a quantum pump 

In this manuscript we address the current noise in an 
adiabatic quantum-dot pump, in which a current is in- 
duced by a slow time-dependent, periodic modulation of 
the system parameters. The expectation value of the 
resulting pumping current is time-dependent and fulfills 
the condition '^^Ia{t) = — e^(n)(t), with the currents, 
Ia{t), into the left and right leads, a = L, R, and the 
charge in the (quantum dot) system, e{n){t) (e < be- 
ing the charge of the electron). When at least two pa- 
rameters of the system are varied in time, charge can be 
pumped through the system in every cycle. In this case, 
as well as when an additional stationary bias is applied, 
the time-averaged current (namely the dc component of 
the current) is, in general, finite. The dc current through 
the system is defined as 

I- f-m- r-liMt)-iK{t)) , (1) 

Jo ^ Jo T ^ 

where r = 2tt/Q is the period and D, the driving fre- 
quency. We are interested in the quasi stationary situa- 
tion with a periodic modulation, where the expectation 
value of the charge on the dot is conserved after one pe- 
riod. Therefore, the time-averaged current / is equal to 
/l — —Ir due to charge conservation. 

The current flowing through the pump is in general 
noisy, due to thermal fluctuations and due to noneqilib- 
rium fluctuations. We are interested in the noise as a 
spectroscopy tool and as a measure of the precision of 
the charge pump. The time-resolved current noise can 
be defined in the symmetrised form as 

Sa^p (t, t') = (<5t (t) Sip it') + sip it') si^ (i)) , (2) 

where the operator Si a (t) = ia {t) — {ia {t)) measures the 
deviation of the current in lead a from its mean value. 
For a time-independent Hamiltonian, this object can be 
written as a function of the time difference, t — t' . For 
systems that are subject to an external time-dependent 
driving it depends in general on two times. We are inter- 
ested in the noise from a long-time measurement, namely 
the zero-frequency contribution to the noise spectrum 

r dt f°° 

Sep = - d{t~t')So.p{t,t') (3) 

Jo J-oo 

=: f-So.p{t) . (4) 
Jo ^ 

Note that for the general noise-power definition, the in- 
tegral in t is going over a large measuring time interval.!^ 
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In the stationary limit, the integrand Sap {t) = Sap (0) is 
time-independent. However, in the case of periodic driv- 
ing, Sap (t) as defined in Eq. Q is periodic in t, and the 
integration interval can be reduced to one driving period 
T. We will in the following calculate the zero-frequency 
noise in timc-dependently driven quantum dots. 

II. MODEL AND FORMALISM 

A. Model Hamiltonian 

As a model system for the study of the pumping noise 
we consider a single-level quantum dot coupled to two 
leads a = L, R, as shown in Fig. [T] The quantum-dot 
system is described by the Anderson model where certain 
parameters in the system are time dependent 

+ Un^ni+ Y ^ (7a (i) 4cQfe<T + h.c.) . (5) 

a— L.R ka 

The first term with the creation (annihilation) operator 
'^akcr i'^aka) dcscribes electrons with spin a and momen- 
tum k in the leads L and R, which can have different, 
possibly time-dependent chemical potentials. The energy 
level on the dot is denoted as et,i(0 — ^(0 "F and 
can be time dependent via driving by a gate. Eventually, 
we consider a finite magnetic field, leading to the Zeeman 
splitting A. The number operator of electrons with spin 
a on the dot is ha- = d^dg., where {d^) creates (annihi- 
lates) an electron with spin a on the dot. Coulomb inter- 
action between electrons on the dot is included within the 
constant interaction model and is here parametrised by 
U. Electrons can tunnel from the leads onto the quan- 
tum dot, where the tunneling matrix element is given 
by the, possibly time-dependent, parameter "fait)- The 
tunnel coupling strength is characterized by the rates 
Ta{t,t') = 27rp„7„(t)7*(t'), with Fait) = Ta{t,t). In 
the following we set /i = 1. 

B. Real-time diagrammatic approach for noise 

In this section we write the zero-frequency current 
noise as well as the current itself based on a kinetic equa- 
tion, which allows for an evaluation in terms of a non- 
equilibrium, real-time diagrammatic approach. In order 
to calculate these expectation values, information on the 
time-dependent density matrix of the system is needed. 
Since we are not interested in the specific time evolution 
of the states of the reservoirs, we trace out the degrees of 
freedom of the reservoirs and calculate the expectation 
values based on the knowledge of the reduced density ma- 
trix of the quantum dot and its time evolution. We are 
then left with a four-dimensional Hilbert space spanned 
by the dot states |x) € {|0), | f), U), |d)}. The state 



associated with an empty dot has energy = 0, the en- 
ergy of the singly occupied dot with spin a =t, i is given 
by E„ — while the energy of the doubly occupied 
dot also includes the Coulomb interaction, Eii — 2e + U. 
We assume that the quantum dot and the reservoirs are 
decoupled at some initial time to; since we are not in- 
terested in transient effects, we will later set this initial 
time to — oo, meaning that it is far away from the mea- 
suring time. The reduced density matrix of the dot at a 
later time t is connected to the density matrix at time 
to via the propagator n(t, tg)- In general, the reduced 
density matrix and the propagator are tensors of rank 2 
and 4, respectively. In the following, however, we con- 
centrate on the time evolution of the diagonal elements 
of the reduced density matrix, summarized in the vector 
Pit) = iPoit),P^it),P^it),Pdit)). This is justified by the 
fact that in the model under consideration the dynam- 
ics of diagonal and off-diagonal elements of the reduced 
density matrix do not couple to each other due to charge 
and spin conservation in tunneling processes. Therefore, 
we write 

Pit)^Ilit,to)Pito) , (6) 

where the propagator Hit, to) is a four by four ma- 
trix connecting the vectors P{t) and P(to)- The (non- 
equilibrium) time evolution can be described using the 
Keldysh contour. It is depicted by an upper (lower) time 
line, representing the forward (backward) part of the con- 
tour, indicated by arrows. Such a contour is shown for the 
calculation of the noise, which will be introduced later, 
in Fig. |3] The propagator fulfills a Dyson equation of the 
form 

n(t,io) = i+ / dti f dt2Witi,t2)n{t2,to) , (7) 

with the kernel W as the self energy, which describes all 
processes coupling the quantum dot and the leads. The 
time evolution of the reduced density matr ix can then be 
written in terms of a kinetic equatiorPsEH 

|p(t)- 1* dt'Wit,t')P(t') , (8) 

where the kernel Wit, t') accounts for transitions between 
different reduced density matrix elements due to tunnel- 
ing. It can be computed diagrammatically; we here use 
a diagrammatic language ,1^^ where the process of an elec- 
tron hopping between dot and reservoir is represented 
as a tunneling vertex on the Keldysh contour. These 
vertices are connected via tunneling lines, which are con- 
tractions stemming from the tracing out of the lead de- 
grees of freedom (due to Wick's theorem, there is always 
a pairwise connection of vertices). In this context, the 
kernel, i.e., the self energy of the Dyson equation 
can be regarded as the sum of all irreducible diagrams, 
see Appendix [C| for more details. 

A similar equation as Eq. ([s]) can be written for the 
time-dependent expectation value of the current, lait) = 
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FIG. 2: Diagrammatic representation of the current expecta- 
tion value {Ia{t)). The current operator originally placed at 
the turning point of the Keldysh contour, can be equivalently 
put on the upper or the lower contour. By summing both 
possibilities, a factor 1/2 is added. 



{Ia{t)), which we consider at a certain measurement time 
t, by introducing the object Wj^. The object Wj^ rep- 
resents the sum of all irreducible diagrams including the 
current operator as an external vertex. In Fig. [2] (left- 
hand side), we draw a graphic representation of W7„. 
The states as shown at the borders of the kernel at times 
t or t' are equal, since, as mentioned above, only diagonal 
elements of P are of interest. The current operator sits 
at the turning point of the Keldysh contour at time t. 
In principle, the Keldysh contour could be continued un- 
til any later time > t without changing the expectation 
value. Here, we choose the shortest representation. The 
aim is to relate Wj^ with the kernel W. We therefore 
cut the contour at the turning point, giving the equiv- 
alent choices to place la either on the upper or lower 



part of the contour (see right-hand side of Fig. [2]) . Since 
both of them are equal, there is a correction factor of 1 /2 
added. Thus we can write 



2' 



dt'Wi^ {t,t')P{t') 



(9) 



with e""" — (1, 1, 1, 1), being a vector representation of the 
trace operator acting on the reduced Hilbert space (in 
Fig. [2] the multiplication of the vector e""" is represented 
as the sum over x')- The current kernel Wi^ as appearing 
in Eq. ^ can be directly derived from the kernel W 
by use of the accordingly modified diagrammatic rules 
shown in Appendix [C| 

We are in the following interested in the zero- frequency 
noise. In order to express the zero-frequency noise based 
on the time evolution of the current operators on the 
Keldysh contour, we concentrate on the object Sapit). 
This is the integrand of the expression for the zero- 
frequency noise, given in Eq. ([s]). We distinguish the 
two cases where the time argument of the current is 
larger than the time argument of the current namely 
t > t', yielding the three integral expressions of the fol- 
lowing formula and the opposite situation t' > t, where 
the current Ip is taken at a later time and which we ob- 
tain by exchanging the indices a and /3, 



Sap (t) =— e lim 

Z to^ — oo 



dtiWi^i^(t,h)P{h)+ dh dt2 dhWi^{tM)^{h,t2)Wi^{t2M)P{h) 

t rt' 



dh dt' I dt3Wi^,{t,t,)P{h)(g,e^Wi^{t\t3)P{t3) 

to J to 



+ (a o /3) . 



(10) 



We discuss in the following the three different integral 
expressions occurring for t > t' , examples of which are 
depicted in Fig. [3] The initial time to before which we 
suppose the system and the leads to be decoupled is not 
taken in the limit tg — > —00 yet. The reason for this is 
that the different integral expressions considered here do 
not converge separately for this limit, since the propaga- 
tor n {ti,t2) does not have a finite support. 



The first term of Eq. ( 10 1 deals with the possibility to 



place both current operators in the same irreducible part 
of the Keldysh contour, captured by the object Wj^j^, 
see Fig. [3] (a). The operator occurs at time t, which is 
the latest time within Sap (t), see Eq. (10), and is placed 



at the turning point of the Keldysh contour (as before, in 
the case of Wi^ see Fig. [2| . The second current operator 
Ip has no fixed position within the object Wj^i^^. We 

here show the example where the current operator Ip is 
placed on the lower contour, however also the diagrams 
where Ip is placed on the upper contour contribute to 



the situation where t > t' . 

In the second term the current operators are in differ- 
ent irreducible parts, which are connected via the prop- 
agator n of the system, see Fig. [s] (b). Also here, we fix 
the operator at the contour turning point at time t, 
whereas Ip can be placed anywhere on the contour (in 
Fig. |3] (b) , the lower contour part was chosen) . 

The last term accounts for the contribution due to the 
product of two separate current expectation values at 
different times, where the tensor product occurs as the 
product of two terms which can be depicted on separate 
Keldysh contours, see Fig.^(c) and (d). The times, t and 
t' , at which the currents, la and Ip, are taken can occur 
in different realizations, where we here consider again the 
case t > t': The time t' can lie within the time interval 
set by the Kernel Wj^ {t,ti), this is shown in the example 
Fig. [3] (c). Note that, furthermore, the relation between 
the times ti and ^3, the times at which the current kernels 
start, is not fixed; in Fig. |3] (c) we show as an example 
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FIG. 3: Diagrams for the cases of having (a) both current 
operators in the same irreducible block, and (b) for having 
them on different blocks. Diagrams (c) and (d) depict the 
last contribution due to the product of current expectation 
values, where for (c) ti < t' and for (d) ti > t' . 



the case where ti > t^. 

Alternatively, the time t' can lie outside the time in- 
terval set by the Kernel Wi^{t,ti), i.e., t' < ti, see 
Fig.[3](d). In this case the third integral expression has a 
similar structure as the second one, depicted in Fig.|3](b). 
We use this fact for the evaluation of the zero-frequency 
noise, extensively discussed in the AppendixjB] where we 
combine these two contributions to one. 

The discussed object, Sa/s (t), is periodic in time and it, 
furthermore, enables us to perform in a subsequent step 
a systematic adiabatic expansion of the noise with re- 
spect to time t. The zero-frequency noise has before been 
studied based on a real-time diagrammatic approach for 
a time-independent system in Ref. 1741 Here, we are in- 
terested in considering the explicit time-dependence aris- 
ing from the pumping cycle affecting the zero-frequency 
noise. The general expression given in Eq. ( 10 ) is the 
starting point of our considerations. The treatment of 



the time dependence is presented in the following sec- 
tions UlC] and EDI 



C. Adiabatic expansion of the current 

We are interested in a situation where the parameters 
of the system depend on time, as indicated in the Hamil- 
tonian in Eq. ([5]). The time dependence is taken to be pe- 
riodic and slow, such that the time scale of the parameter 
modulation is much larger than the lifetime of the states 
in the system. We aim to write Eq. ( 10 1 for the limit of 



adiabatic driving. We can, hence, perform an adiabatic 
expansion of the kinetic equations introduced above, fol- 
lowing the lines of Ref. _27, This expansion contains in 
lowest (zeroth) order in the frequency the " instantaneous 
contribution" of all terms contributing in Eq. ( 10 1 which 



arises from a situation where the system can be taken 
to always follow the parameter modulation. This means 
that all parameters are frozen to time t. An "adiabatic 
correction", in first order in $7, accounts for the actual 
delay with respect to the instantaneous solution. In the 
following we indicate the contribution in zeroth order of 
the frequency expansion with the superscript (i) and we 
denote the first-order correction in frequency by (a). 

In this spirit, the kinetic equation for the adiabatic 
expansion of the reduced density matrix of the dot. 



Pit) ^ Pt 
mated as 



Pt°'\ starting from Eq. 



d 
dt 



= {WP}[''> 



IS approxi- 

(11) 
(12) 



The instantaneous contribution and the adiabatic correc- 
tion to the reduced density matrix are found as solutions 
of these two sets of equations, which contain the instan- 
taneous contribution and the adiabatic correction to the 
kernel, W W'^'^ + W'-''\ as well. We use a short no- 
tation with brackets, which is generally defined for two 
arbitrary objects A and B (be it any of the kernels, the 
density matrix, or also products of several objects) as 



(a) r(*) 



A^'B. 



dAY'B. 



(13) 
(14) 



This expression can be readily derived from a general ex- 
pression for the frequency expansion of a convolution, as 
previously done by Kashuba et al. Ref. |80l see also Ap- 
pendix [X] The Laplace transforms of the instantaneous 
kernel and its adiabatic correction, Wi^^'^\z) = d{t — 
i')e-(*-*')l^/^/'^)(< - t') = /^^di'e"-(*-*')M^/*/"'(t - 
t') occur in these equations, where the short notation 
without the 2-argument denotes the limit VFj 



lim^ 



(i/a) 



(z) and dW^ 



lim^ 



[dW,^'\z)/dz]. 



The kernel W; ' can be constructed in a straightfor- 
ward manner by taking the kernel of the stationary sys- 
tem and inserting the time dependence parametrically. 
This is indicated by the subscript t. The adiabatic cor- 
rection to the kernel contains the effect of the time- 
dependent parameters changing in time along the contour 
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in first order in Q and is obtained according to the rules 
given in Appendix [C] see also Ref. [27l The finite width 
of the kernel enters this equation through the derivative 

In order to justify the adiabatic expansion of the mas- 
ter equation, one has to take into account two relevant 
time scales. For one there is the dynamics of the density 
matrix, dominated by the tunnel coupling to the reser- 
voirs, scaling with the order of magnitude of W itself. 
And there is the finite support in time oiW{t — t'), which 
needs to be sufficiently small to perform the expansion 
in 17. This second quantity is represented by the width 
of the kernel dW, and corresponds to the dynamics of 
the reservoir. In this spirit, the adiabatic expansion is 
justified if both the tunnel and the reservoir dynamics 
are much faster than the external driving. 

Due to the slow parameter modulation, also the cur- 
rent through the system has an instantaneous contribu- 
tion and a first-order-in-frequency, adiabatic correction. 
We will refer to this first-order correction as the pump- 
ing current. In analogy to the kinetic equation one can 
compute the current of the system by the respective adi- 
abatic expansion of the current kernel and the reduced 
density matrix 



(a) 



(15) 
(16) 



The current kernel contributions W^j*^"' and dWj'' ^ are 
constructed according to the rules in Appendix[C] The re- 
duced density matrix, obtained from Eqs. ( [TT| ) and (12 1, 
enters the formula for the current expectation value. Im- 
portantly, if no bias voltage is applied, the instantaneous 
current is always zero. In this case the pumping current, 
/q"\ becomes the dominant contribution. We are in the 
following interested in the situation of zero bias, namely 
the regime of pure pumping as well as in the situation 
where a pumping current flows on top of an instanta- 
neous current due to a non-equilibrium bias. 



D. Adiabatic expansion of the noise 

The main purpose of this paper is the calculation of the 
zero-frequency current noise of an adiabatic quantum- 
dot pump. The aim of this section is to present the 
zeroth-order and first-order contribution in the driving 



frequency to the zero- frequency noise, s'^l and S^g, go- 



a/3 



ing along with the pumping current presented in the pre- 
vious section. For this purpose we expand the elements 
of the periodic function Sap (t), given in Eq. ( 10 ), around 
the reference time t. In a first step we consider the instan- 
taneous contribution to the different terms, only. The de- 
tailed derivation is given in Appendix [B] We obtain for 
the zeroth-order contribution in the driving frequency to 



this function 



s': it) 



{Llpf^^^ + {a ^ P) , 



(17) 

where we describe the different elements in the following. 
The first term contains the kernel W/^/^ involving two 
current operators; the functional form of this expression 
is analog to the one found for the instantaneous current. 
The second term, containing current kernels from differ- 
ent parts of the contour is a combination of the second 
and part of the third expression of Eq. (10 1. Note that 



an objec t related to the propagator connects the two ker- 
nelPEEl 



n(t,i') ^u{t,t') - p{t) 



(18) 



This object has a finite support and is, hence, zero for 
large time differences, t — t' ^ oo. Its further proper- 
ties and its evaluation starting from a Dyson equation 
are discussed in Appendix |d] The important trait of 11 
in order to justify the adiabatic expansion, is that the 
width of n needs to be sufficiently small. The relevant 
time scale determining the support of 11 is given by the 
dynamics of the density matrix, and hence, scales with 
the time scale determined by the inverse of the tunnel- 
ing rate. Therefore, there is no new relevant time scale 
appearing with respect to the ones discussed before, and 
the same conditions that justify an adiabatic treatment of 
the kinetic equation, Eqs. (15) and (16), justify likewise 



the expansion of the noise for small pumping frequencies 

The condition needed to evaluate the instantaneous 
contribution to the object 11 in Laplace space (again, 
taking the z — > 0+ limit) is 



[Wli] 



{dWP} 



ii) 



(19) 



with the extra condition e'^II 



0. 



Finally the third term of Eq. ( 17) contains the remain 



ing part of the product of two current expectation values 
at different times. In this third integral expression, the 
object 



t2 



dhWiStM)P{h) , (20) 



occurs, which is an object depending on two time- 
arguments and differs from the definition of the current, 
Eq. ([9]), by the integration interval. The Laplace trans- 
form of this object is discussed in detail in Appendix [B| 
Its zeroth-order contribution in the driving frequency 
for z -> 0+ is given by Yiia^^f)+ ia\t, z) = —dla\t), 
where we introduced the abbreviation 9/q (t) = 
le^dwi%P,^'\ 
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Integrating the function 5*^'^ of 



Eq. ^ 



over one pe- 
riod in order to obtain the zero- frequency noise as de- 
fined in Eq. ([s]) dehvers a time-averaged version of the 
stationary zero-frequency noise considered by Ref. I74I81I 
We will show later, that if the time-dependent driving 
of the system occurs with large amplitudes, already the 
noise in zeroth order in the pumping frequency can devi- 
ate strongly from the one of the corresponding stationary 
system. 

The important object for this paper is the pumping 
noise, namely the first-order in the pumping frequency 
correction to the noise, arising from the slow driving. 
Only if there is a working pump, i.e., only if there are 
two time- dependent parameters enclosing a finite area 
in parameter space in one period the pumping noise is 
non-zero. With the help of the notation using braces, 
introduced in the context of the kinetic equation in the 
previous section, the elements of the pumping noise can 
be written as 



{a) 



(a) 



-2 



(21) 

The first two elements can be readily identified as the 
adiabatic corrections to the first two terms in Eq. (17 1. 
In addition to the adiabatic correction to the reduced 
density matrix and the current, they contain the adia- 
batic correction and the first derivative with respect to 
the Laplace variable of the object 11 related to the prop- 



agator. The recipes to evaluate dll^ 
by 



J Tf('') 

and 11( are given 



d{WU}f' = n 



{wn} 



t 



P. 



(a) 



1 



2 [d'WP]^ 



-{dwp} 



(a) 



with the normalization conditions e'^'ll 



(a) 



(22) 

e^ . 

(23) 

and 



In these expressions we find not only the first but the sec- 
ond derivative of the current kernels with respect to the 
Laplace variable z, d'^w'f \ — \mi^^^+ dlw[^\{z). With 
these sets of equations the current and the zero-frequency 
noise can be evaluated as soon as explicit expressions for 
the different kernels are found as explained in the follow- 
ing, see Appendix [O for detailed rules. 



E. Perturbation expansion in the tunnel coupling 

In the following we assume that the coupling between 
the quantum dot and the reservoirs is weak. This allows 
us to perform a perturbation expansion in the tunnel 
coupling F on top of the adiabatic expansion performed 
above. This perturbation approach is a useful approx- 
imation as long as the condition F <C fceT is fulfilled, 
where T is the temperature and fee the Boltzmann con- 
stant. Within this perturbation expansion, the Coulomb 
interaction and the non equilibrium due to an applied 
bias can be treated without any further approximations. 
The kernels describing tunneling events between quan- 
tum dot and reservoirs are always expressions in at least 
first order in the tunneling coupling strength F. We write 
their expansion in general as 



^ * X.t 



(27) 



where the index x stands for the different kinds of ker- 
nels, such as W , the current kernel Wi^ , or Wi^ii^ . While 
expressions containing the kernel only in the lowest-order 
expansion in the tunnel coupling are commonly referred 
to as the sequential tunneling contribution, we here go 
beyond this limit by taking into account quantum fluc- 
tuation effects put to second order in the tunnel coupling 
F. 

In order to find the instantaneous and adiabatic con- 
tributions to the reduced density matrix, the current and 
the zero-frequency noise in lowest and in next-to-lowest 
order in the tunnel- coupling strength F a rigorous or- 
der by order expansion of the kinetic equations given in 



Eqs. (11) & (121, the currents from Eqs. (15l & (16), the 



The adiabatic correction to the third term of Eq. (17) zero-frequency noise from Eqs. (17 1 & (21), as well as 



contains again the object /. In order to evaluate the adi- 
abatic correction of the bracket the knowledge 

of the derivative with respect to the Laplace derivative 
of is needed 



(24) 



This is discussed in further detail in Appendix |Bj Two 
more definitions are needed to evaluate the contributions 
to the adiabatic correction to the bracket containing /, 
which are 



2r(») 



it) =36^5^0/'^ 



it) =^e'' {dWj^P} 



(a) 
t 



(25) 
(26) 



the propagator functions entering the noise expressions 
as shown in Eqs. (18), (22) and (23) needs to be per- 
formed. The results are discussed in the remaind er of 
this paper. It is known from previous calculations j ^^ * ^^ 
that the lowest-order contribution to the reduced density 
matrix in the instantaneous limit is the zeroth order in 
F, denoted by p(*'"), while the adiabatic correction starts 
in minus first order, p(°^~i). This is consistent with the 
fact that fi/F ^ 1 in the adiabatic limit. The domi- 
nant contribution to the instantaneous current, l'''''^\t), 
is a first-order in F contribution, while the pumping cur- 
rent starts in zeroth order in F but is proportional to 
the frequency, /'"'"^(f). Likewise, the instantaneous zero- 
frequency current noise begins in first order in F, sj^'^\ 
where the pumping noise contribution starts one order 
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lower, S, 



(a,0) 

a/3 



In the remainder of this paper we apply 



these general formulas to the single-level quantum dot, 
introduced in Secjllj for a slow, time-dependent driving, 
in the zero- as well as in the finite-bias regime. 



III. PUMPING NOISE 

In this section we want to study the zero-frequency 
noise of a single-orbital quantum dot, as introduced in 
Sec. |II A[ We consider the case of zero bias voltage, where 
a finite pumped charge is due to the modulation of at 
least two time-dependent parameters. In the absence of a 
bias the instantaneous current is zero for all times. How- 
ever, a finite charge can be pumped through the dot, i.e., 

the dc pumping current is non-zero, if there is a fi- 
nite phase difference between the two parameters. These 
parameters can be the level position of the dot and the 
tunnel-coupling strengths to the different contacts, which 
can all be driven by externally applied gates. In the fol- 
lowing we write these parameters as X{t) = X + SX{t), 
where X denotes the time average of X{t). 



A. Instantaneous contribution to the noise 

While the instantaneous current always vanishes at 
zero bias, thermal noise is present, even without the mod- 
ulation of external gates. We start the study of the noise 
in the unbiased case by this instantaneous contribution 
and show to which extent this thermal noise reveals fea- 
tures of the time-dependent driving. In lowest order in 
the tunnel coupling, we find 



'lr 



_4e2 f ^rLlR^^^.^I^O) 



(28) 



Note that we can express the noise exclusively in 
quantities, which are local in time and which charac- 
terize the properties of the quantum dot.f^ Namely, 
Eq. (28) contains the relaxation rate of the charge 



on the dot in lowest orde r in the tunneling, A^^' — 
r[l + /(e) - /(e-H [/)], 1221841 with the Fermi function 
f{E) = l/[l+exp{/3E)] and /3 = l/fcer, as weU as the in- 
stantaneous charge variance on the quantum dot, defined 
as 



(Ac)!-") 



2\(^0) 



/ 2\ I 



(c)J 



(29) 



with the definition (c) (n) —1, where (n) is the average 
charge on the dot. The quantity (c) represents an aver- 
age electron-hole occupation, namely it is 1 if the dot is 
doubly occupied, and —1 for an empty dot, i.e., there are 
two holes. Otherwise, i.e., for one electron on the dot, c 
is zero. We choose this representation since the quantity 
(c) appears in noise expressions much more often than 
(n) alone. The reason for this is that the noise is related 
to charge fiuctuations. We will see in the following, that 



the noise reveals interesting symmetries with respect to 
the electron- hole symmetric point. The explicit result for 
the quantum-dot charge variance can be given as 



= 1- 



A 



2r 



1 - 



(30) 



with (c)f = [1 - /(6) - fie + U)] I [1 + /(6) - /(6 + [/)]. 
The result for the instantaneous zero-frequency noise 
is shown in Fig. [4] (a), for the case of time-dependent 
energy level and left tunneling rate, t — Z be sin(ilt) 
and Fl = Fl -f 5Fl sin(J7i -I- (^), where pumping am- 
plitudes are small, i5FL,5e ^ f. The instantaneous 
cross-correlation is always negative, as required.!^ For 
small amplitudes it shows a resonant contribution at 
e « 0, when the addition energy of adding an electron to 
the empty dot is equal to the electrochemical potential, 
and at e ss — C/, when a second electron can be added. 
These two contributions coincide, when the Coulomb 
interaction is zero (see blue, dashed line). One can 
observe, that the single peak for J7 = is not the sum 
of the two peaks for finite J7, but is reduced by a factor 
of roughly 0.7. In Eq. (28 1 we see that this arises from 



the fact that the tunneling rate prefactor is modified by 
the Coulomb interaction through Ai^''/F. For ?7 7^ 0, 

there is an increase of the rates, since IV > 1 when 
e is within the range from —U to 0. This leads to an 
effective increase of the total current fluctuations for 
flnite U. 

Since no bias is applied, the instantaneous part of the 
noise is a pure thermal noise. Equation (28) can be un- 



derstood as a time-averaged fluctuation dissipation the- 
orem for it can be related to the instantaneous linear 
conductance as 



c(i) _ 



-4.kuT 



dt 



G 



where the linear conductance is defined as 



G 



dV 



(31) 



(32) 



This result holds also in next order in the tunnel cou- 
plin g, wit h corrections to the linear conductance of the 
fornPMl 



de 



(33) 



with the lowest-order conductance G^^\ which can be 



directly extracted from Eq. (28), and the definitions 



e-f [/)-</>(£) 
ar= 1 |:a-<^'(e)-0'(e + C/) 



(34) 

l-f{e)-f{e + U) 

1 - /(e) + ./(e + U) ■ 
(35) 



We defined 0(w) = ^Re^- -h if^) and 0'(w) = 
where ^' is the digamma function. The term 
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0.002 L 



FIG. 4: Plot for (a) the instantaneous and (b) the adiabatic 
zero-frequency noise as a function of the mean energy level e. 
The graphs are drawn for [7 = 0, 20^3^ (blue dashed, black 
solid). The left and right tunnel coupling rates are Fl = 
0.3F and Fr = 0.7F. The other parameters are F = Q.^k^T, 
the pumping amplitudes (5e, (5Fl = O.IF, the phase difference 
between the pumping parameters = 7r/2, and Q. = O.IF. 




10 ksT 



FIG. 5: The instantaneous noise (black solid) and the sum of 
instantaneous and adiabatic (red dashed) as a function of the 
mean energy level e. The splitting of the resonances is visible 
for 5e > ksT, here 5e = AksT and ST = 0.3F. The other 
parameters are U = 2Qk^T, Q. = F/4, and F = O.SfeBT. 



a denotes the renormalization of the energy level due to 
quantum fluctuations, e — > e + a. The function uy on 
the other hand is related to the r-renormalization. The 
prefactor depending on the Fermi functions relates to the 
fact that (Tr contains the weighted F-renormalizations at 
different resonances, occurring with opposite signs. For 
details see Ref. |86l Also, see Refs. |85| and [86| for the 
remaining analytic expressions for the cotunneling con- 
tributions to the conductance. 



Importantly, it is the time- averaged conductance which 
determines the instantaneous zero-frequency noise. This 
becomes clearly evident in the black (full) line of Fig. [S] 
for the case of large amplitudes 5t ^ k-oT, where the 
time dependence of the level position can be seen from 
the splitting of the resonances. This feature is due to 
the fact that when the amplitude of the level energy 5e 
is much larger than ^bT, the quantum dot can be at 
resonance for a considerably large time, even though the 
mean energy level e is far from resonance. When the 
distance from e to the resonance is equal to the amplitude 
at harmonic driving, the dot level is at resonance at the 
turning point of the sinusoid, hence the dot stays in the 
resonant regime for a much longer time than if the mean 
energy were exactly on resonance. 

Moreover, for large amplitudes, the instantaneous ther- 
mal noise part gets decreased considerably, as can be 
seen when comparing Fig. |4][a) with Fig. [5j due to the 
fact that the addition energies are not always close to 
the resonanc es dur ing the whole pumping cycle. We will 
show in Sec 



HID 



that 5^]^ can even vanish for certain 



pumping schemes with large modulation amplitudes. 



B. Adiabatic pumping with e{t) and Vc,{t) 

While an instantaneous contribution to the noise is 
also present in a time-independent system, an adiabatic 
contribution is uniquely due to a finite time-dependent 
modulation. Therefore the characteristics of the pumping 
noise contain more information about the actual pump- 
ing process than the charge only, as we will show in the 
following. 

By assuming zero bias, but an arbitrary combination of 
time-dependent parameters, we find for the leading order 
in the tunnel coupling for the first-order in O correction 
to the noise (i.e., the pumping noise) 



7 r2 



(36) 



Importantly, also for the case of time-dependent driving, 
it is possible to express the zero-frequency noise exclu- 
sively in terms of observables in the dot which are local 
in time, and their derivatives. We immediately see that, 
similar to the pumping current,'^ the pumping noise in 
the sequential tunneling regime can be non-zero only if e 
is one of the time-dependent parameters, such that the 
occupation number {c)[^'^^ changes in time, too. For the 
remainder of Sec. |IIIB| we choose Fl as the second time- 
dependent parameter. 



Comparing Eq. ( 28 ) and Eq. ( 36 ) , we see clear similar- 



ities between the structure of the pumping noise and the 
instantaneous noise, which exclusively stems from ther- 
mal fluctuations. However, the adiabatic noise is due to 
a working pump (i.e., two independent driving parame- 
ters) and is, therefore, a non-equilibrium effect. Hence 
unlike in the instantaneous case, where the system is as- 
sumed to be in equilibrium for each point in time, one 
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cannot, in general, expect the fluctuation-dissipation re- 
lation, Eq. (31 ), to hold for the pumping noise. However, 



within the framework of a scattering-matrix approach 
(valid for U = 0) we are able to show that the equi- 
librium relation can indeed be extended to first order in 
n in the following way 



^ dt 



g|°^(c/ = o) 



(37) 



Note that this relation holds for arbitrary orders in F, 
provided that the temperature is still larger than the 
driving frequency fl. Equation (37) is a consequence of 



the special case where the system both is non-interacting 
and adiabatically driven (i.e., close to equilibrium). Re- 
garding the formulation of general fluctuation relations 
out of equilibrium, we refer to Refs. l50|87ff9Tl 

In the following we show that in the presence of 
Coulomb interaction, a correction term appears, and the 
relation (37) breaks down already in first order in F. For 
the interacting system, the pumping noise, Eq. ( 36 ) , can 
be expressed as 



,(a,0) 

'lr 



_4fcBr /'"^g("^")(o 







" dt FlFr 
T F2 



AW-f) (c)^°) 



d_ 
di^ 



Am 



(38) 



It can thus be separated into a part proportional to 
C(i,o)^ which represents an adiabatic correction of the 
equilibrium fluctuation dissipation theorem, and a cor- 
rection term. Importantly, for a non-interacting system, 
xi^^ = F, and therefore the correction term vanishes. 
Hence it is the Coulomb interaction that prevents us to 
write down all the noise corrections in terms of the con- 
ductance. This is a strong indication that 5'^'^^'' is not 
uniquely of thermal origin for an interacting system. 

The adiabatic noise is shown in Fig. [4] (b) for the case 
of small pumping amplitudes. We see for one that the in- 
stantaneous noise is much larger than the adiabatic one. 
However, an important distinguishing feature is their 
symmetry: whereas the instantaneous cross-correlation 
is always negative, its flrst-order in correction can be 
both positive or negative, going along with a sign change 
at each resonance. Therefore, the pumping noise can 
lead to an enhancement as well as a reduction of the to- 
tal noise. In addition, the pumping noise has a node 
at the electron-hole symmetric point I — —U/2 and is 

antisymmetric arround this point (in contrast to S^^, 
which is symmetric). For small amplitudes, as shown in 
Fig.|4](b), this leads to a shift of the purely instantaneous 
noise contribution, which takes place in the same direc- 
tion for both the contribution at e and at e + C/ w 0. 
This shift direction of both contributions depends further 
on the phase-difference of the driving parameters and the 
coupling asymmetry of the tunnel barriers.'sS 

Considering the adiabatic correction, in the regime of 
large amplitudes it results in altering the heights of the 
two parts of the separated peaks (as a reminder, the 



Fl 


-Fr = 


-o.4r 


Fl 


-Fr = 


o.2r 




FIG. 6: Plot of the adiabatic Fano factor i^'"' as a function of 
the mean energy level ?. The two curves are for (Fl — Fr) /F = 
—0.4,0.2 (black, red). The other parameters are F — O.bksT, 
the pumping amplitudes 5e, 5Fl ~ O.IF, the phase difference 
between the pumping parameters ip — n/2, and Q — O.IF. 



peaks separate due to the large amplitudes), rather than 
in a shift. Again, which of the contributions gets in- 
creased or decreased depends on the coupling asymmetry 
and the direction of the pumping. 

In order to characterize the zero-frequency noise due 
to pumping with respect to the average current pumped 
through the system per period, /'■°'', we define an adia- 
batic Fano factor 



(a) 
LR 



(39) 



Since there is no instantaneous current in the absence 
of a bias voltage, and thus / =0, this quantity can 
be interpreted as an adiabatic correction to the ordinary 
Fano factor. It gives information about, e.g., whether 
there is a time-averaged pumping current that is free of 
pumping noise (i^(°) is zero), or there is pumping noise 
in absence of a pumped charge (F^"^ diverges). 

Note that the direction of the dc pumping current and 
the sign of the related pumping noise are highly sensitive 
to all system parameters, which is in contrast to station- 
ary transport. Therefore our interest in the adiabatic 
Fano factor F^"^ is to examine the relation of the cur- 
rent and the noise including also their respective sign, 
differing from the motivation for the original Fano factor 
namely to study the noise with respect to its Schottky 
limit. 

Importantly, F^"-* is independent of the pumping di- 
rection, because both /^°^ and change sign for a re- 
versed pumping cycle. This means that the sign of F^^^ 
can be related to intrinsic characteristics of the system 
(as we will see in the following). The adiabatic Fano 
factor, F^"-\ is shown as a function of the average level 
position in Fig. [6) We focus on the case of small ampli- 
tudes of the parameter modulation, where we can derive 
meaningful analytic results. In a bilinear expansion of 
the pumping amplitudes, we find for the adiabatic noise 



(a,0) 
LR 



^^,M Fr(Fl-Fr) d 
T de 



(Ac) 



(40) 
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with the area SA = JJ" dtTi^e, which is enclosed in pa- 
rameter space in one pumping cycle. As for the time- 
dependent system quantities, the bar over an arbitrary 

fmiction, e.g. (Ac)[ ' \ always indicates that only av- 
erage parameters appear. For Fl = Fr, we find that 
S'l'^°\ in contrast to the pumped charge, is zero. Note 
that while the instantaneous noise and the pumped 
charge have even parity with respect to the particle- 
hole symmetric point, the pumping noise shows odd par- 
ity, 4r°' (-£ - U/2) = -S'^"/' (e + C//2), leading to the 
shift of the total zero-frequency noise, discussed before. 
In lowest order in the tunnel coupling and for small am- 
plitudes, Se ^ ksT and ST ^ F, the adiabatic Fano 
factor can be approximated as 



is due to the level renormalization term, only. 



" Fl - Fr d 



dt 



y (i,rcn) 



F 



The renormalization of (c) is (c)^*''^™'' — a-^[c)__ 
Here, we present the contribution to the zero-frequency 
noise in second order in the tunnel coupling. The adia- 
batic correction in second order is given as 



(42) 



dt FlFr d (j,ron) 



This splitting up of terms is motivated by Eq. (42). Anal 



ogously to the pumped charge, the first contribution is 
related to the level renormalisation, here specifically, the 



Ft,-F 



R d(Pey 



(c) 



(i,0) 



renormalized charge variance (Ac 



(2,rcn) 



R 



d 



(41) 



However, there is an additional correction term 



dm 



'-((a,corr) 

'lr 



dt FlFr 1 d 



Evidently, the adiabatic Fano factor scales directly with 
the coupling asymmetry. This can be seen in Fig. [6j 
where we show the Fano factor as a function of the av- 
erage energy level e for different (i.e., opposite) coupling 
asymmetry. The latter factor can be understood as fol- 
lows. If the dot level is far from resonance and the charge 




(»,o) 



Od 



is exponentially suppressed, (c) 



„±l3e where we used a' = -^a as well as S 



the adiabatic Fano factor is constant with respect to the 
mean level energy F*^"^ = =F (Fl — Fr) /F. If the quan- 
tum dot is in resonance, i.e., the charge depends lin- 
early on the mean level energy, (c)^ ' ~ ±e, we can 
observe a node in F^'^K Furthermore, the Fano factor 
changes sign at the electron-hole symmetric point, since 
the pumped charge and the pumped noise have different 
symmetries with respect to the latter. Hence, the off- 
resonant plateau heights are directly equal to the cou- 
pling asymmetry, whereas the nodes occur at the res- 
onances, and the electron-hole symmetric point. Thus, 
F'"-' displays a distinct feature for the electron-hole sym- 
metric point, which cannot be observed by the pumping 
noise directly, since the noise signal itself is exponentially 
suppressed in this regime. The nodes in F^"'^ correspond 
directly to the nodes of 5'^^ . If F^°-'> = then there is 
no pumping noise even though there is a finite pumped 
charge. In the case considered here, there is no possibility 
to realize the opposite situation of having finite pumping 
noise in absence of pumped charge, in which case the adi- 
abatic Fano factor would diverge. We will encounter this 
case of finite pumping noise with zero pumped charge in 
Sec.HVl 

As a next step, we study the adiabatic noise correc- 
tion in next higher order in F, i.e., in the regime where 
cotunneling contributions as well as interaction-induced 
renormalization effects to the bare dot parameters arise 
due to quantum fluctuations. It has been shown^'' that 
the average pumping current beyond sequential tunneling 



(44) 



i-f{<^)-f{<^+u) 

l-f{e)+f{e+U)- 

The function S has odd parity with respect to the 
electron- hole symmetric point, and gives evidence of 
either the electron-like or hole-like nature of the dot 
spectrum.!^ The correction can be expressed in terms 
of cotunneling transition rates, for one there are rates 
related to transition where initial and final states are 
the same w = PFo-fcr-fO = Wa-^o^^, w = W^^d^cr = 
Wd-^CT-i-d, and also transitions between zero and double 
occupancy occur, Wgd and Wdo- Their explicit form can 
be found in Ref. HI For [/ = 0, we find S^^^^^ = 0, 



in analogy to the dc pumping current / that van- 
ishes also. Moreover, when consulting Eq. (44) we see 

that this correction term is proportional to and 
therefore in spite of the appearance of cotunneling terms, 
the adiabatic noise is still a purely resonant feature, i.e., 
it is exponentially suppressed far from resonance, quite 
like the pumped charge. Moreover one can show starting 
from Eq. ( 43 1 , that in the case of symmetric coupling, 



Fl = Fr, all adiabatic terms but S'[r™"^ vanish. The 
behavior of this correction term is therefore quite differ- 
ent from the other contributions of the pumping noise 
which all vanish for symmetric coupling. We consider 
again the adiabatic corrections of the Fano factor, see 
Eq. ( 39 ) . For a sufficiently asymmetric tunnel coupling, 



the only visible feature is the shift due to the level renor- 
malization, since the correction is much smaller. If we 
however approach the symmetric coupling case, we find 
that F*^°^ does not go to zero, and the only visible fea- 



ture comes from 5, 



(a, I 
LR 



Importantly, S'lr'^"'^'^'' is even 
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the leading contribution to the pumping noise, since also 
the lowest-order contribution vanishes for symmetric cou- 
plmg, SIj^ '. 

Note finally, that the adiabatic correction of the Fano 
factor gives no account for the quality of the pump. For 
this purpose, one needs to include the contribution from 
the instantaneous, thermal noise. When choosing only e 
and Fl as pumping parameters, this is the dominant con- 
tribution, since pumping takes place close to resonance. 
However, when including Fr as a third time-dependent 
parameter, it is possible to reach the quantized pumping 
regime. 



This case we will address in Sec. Ill D 



C. Pumping with the barriers only 

As shown in the previous section, the pumping current 
and the pumping noise in lowest order in the tunneling 
coupling vanish if the level position is not time depen- 
dent, e{t) = Z. The first non- vanishing contribution of 
time-averaged pumping current and pumping noise for 
pumping with the barriers only, is due to second-order 
processes in the tunnel coupling. We use Eqs. (43 1 and 
(44), which are valid also for a modulation of FlJF) and 



rR(t) to evaluate the pumping noise now for pumping 
exclusively via the barriers. Interestingly the correction 
term, Eq. (44), vanishes in this case since it requires 



rent, when pumping with the barriers only, it has been 
found that pumping is uniquely due to level renormaliza- 
tion effects. For the only remaining contribution to the 
zero-frequency pumping noise, we find 



'lr 



2e' 



dt FlFr d 
'V di 



(Ac) 



(45) 



Note that in this situation, the pumping noise is one order 
of magnitude smaller than when pumping also with the 
level position. However, as we have pointed out already, 
the lowest-order contribution of the pumped charge is 
also zero when e is constant, hence the adiabatic Fano 
factor is still of similar order of magnitude as when pump- 
ing with e and Fl. Furthermore, one can show that when 
focussing on weak-amplitude pumping and performing an 
expansion in the pumping amplitudes i n b ilinear order, 
F'") is equal to the expression in Eq. (41). Therefore, 



in the present case of weak pumping and no bias, the 
shape of F^°-'> is robust with respect to the choice of any 
two parameters out of e, Fl and Fr. Importantly, the 
Fano factor contains only the intrinsic properties of the 
system, and the pumping scheme does not obscure F^°'\ 
Keep in mind that this finding is true if the pumping 
consists of only two parameters that drive the system. 



D. Quantized pumping 

In the previous sections, we discussed predominantly 
the case of weak pumping, or when discussing strong 



pumping, we restricted ourselves to only two time- 
dependent parameters, with which there cannot be quan- 
tized pumping. Here, we want consider exactly this 
regime where we transport one electron per cycle. This 
regime is of interest for a quantum standard for the cur- 
rent, because in the regime of qu antized pumping, the 
zero-frequency noise vanishes.lSHZl 

In order to reach a quantized pumping regime for our 
model of a single-level quantum dot, all three parame- 
ters e, Fl, and Fr are required to be time dependent, 
which is in contrast to, e.g., a double quantum-dot sys- 
tem, where two time-dependent parameters are sufficient 
to obtain quantized pumping.EES xhe quantized single- 
quantum dot pump on the other hand is most efficient 
when there is a phase difference of each tt/2 between Fl 
and e as well as between e and Fr , with a total phase shift 
of TT between Fl and Fr . The only possibility to reach 
one pumped charge per cycle at the resonances is when 
the modulation of the F's is such that each F^ closes com- 
pletely at one point in the cycle. Thus, one achieves that 
the pump actually receives an electron from one lead, and 
half a cycle later, it can be released in the other. The fact 
that there is a phase shift of tt between the two F's also 
assures the validity of the adiabatic condition, since the 
sum Fl -I- Fr remains constant. 

We show the result of the number of pumped charges 

r/l"'"^ in Fig. and the corresponding total noise 



7^ 0. Indeed, already for the average pumped cur- 5lr 



MA) 
^LR 



5, 



(a,0) 
LR 



in Fi 



ig- Rb)- 



We see that with 



increasing amplitude of the leveTenergy Se, we reach the 
regime of quantized pumping, / = ef2/27r. When one 
reaches either the resonance at the addition energy of 
one electron on the dot, e = 0, or the one to add a sec- 
ond electron, e = —U, there is one charge transferred 
per cycle, for sufficiently large level energy amplitude. 
Considering the total noise ^lr in Fig. [7](b) we indeed 
find that the noise vanishes in the regime of quantized 
pumping. 

At edges, i.e., when the average energy level e is ±6e 
away from the resonances, there is a remaining thermal 
noise contribution. Here, the dot level is brought close to 
resonance right at the turning point of the sinus modu- 
lation of e where both tunneling barriers are open. Con- 
sequently, there are non-vanishing thermal fluctuations, 
and there is no perfect transmission of exactly one charge 
per cycle. 



IV. NOISE INCLUDING FINITE BIAS 
VOLTAGE AND MAGNETIC FIELD 

We considered until now pumping at zero bias volt- 
age. However, also the study of the current induced 
by a time-dependent bias opens new avenues for trans- 
port spectroscopy. In this section, we therefore generalize 
our previous results by allowing for a finite voltage bias, 
eV{t) = /XL(t) — AtR(i), with the aim to consider adiabatic 
pumping by the modulation of the gate and bias voltage, 
e{t) = € + (5esin(m) and eV{t) = eV + e6Vsm{nt + ip). 
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FIG. 7: The number of pumped charges r/j"''''' (a) and the 
total noise cross-correlations Slr as a function of the average 
level energy e. The time-dependent parameters are chosen 
as e(t) = e + 5es\n(Q,t), Tl = r/2 [1 + sin(r2t - 7r/2)], and 
Fr = r/2[l + sin(m-|-7r/2)]. The amplitude of the level 
modulation is 5e = {SF, 25r, 75r} for the black (solid), blue 
(dotted), and red (dashed) line respectively. The other pa- 
rameters are U = lOOfcaT and k^T = 2r. 



This specific pumping prescription has a very prom ising 
application as a spectroscopy tool. It has been showrpSlSSI 
that when pumping with bias and gate voltages, the 
adiabatically pumped charge is uniquely due to a finite 
Coulomb interaction. If the time average of the bias is 
finite, the occurrence of a pumped charge reflects the 
sign of the coupling asymmetry to the leads, Fl 7^ Fr; 
at zero average bias, the pumped charge contains infor- 
mation about the degeneracy of the ground state. This 
latter point could be proven by a study of the pumped 
charge in the presence of a magnetic field, which lifts the 
spin degeneracy of the system. 

The pumping noise turns out to have a different be- 
havior, containing information on both the pumping and 
the instantaneous current. In particular, we can show 
that a pumping noise can persist in the absence of an 
adiabatically pumped charge. In this case the pumping 
noise indicates that the pumping current cancels in aver- 
age but not at every instant of time. In order to study the 
full parameter space needed for the type of spectroscopy 
introduced in Refs. 28 29j we allow in the following for a 
finite external magnetic field as well. After deriving the 



general formulas for the pumping noise in Section IV A 



we consider the specific case of time-dependent bias and 
gate voltage, for which the pumped charge reveals fea- 
tures due to Coulomb interaction in Section HVBI 



A. General expressions for current and noise 

We start by adding a finite bias V as well as a Zeeman 
splitting A due to a finite magnetic field B to the system 
and study the pumping current and the instantaneous 
and adiabatic zero-frequency noise for a general set of 
time-dependent parameters e, Fq,, A or V . As a conse- 
quence of the bias, there is also an instantaneous current 
flowing in addition to the pumping current, due 
to time-dependent driving. Likewise, on top of thermal 
and pumping noise, shot noise emerges in the large-bias 
regime. 

The magnetic fleld resulting in a Zeeman splitting of 
the energy level adds additional structure to the system 
and, together with the finite Coulomb interaction, ef- 
fects a coupling of the charge and spin dynamics, as will 
be discussed in the following. This coupling of charge 
and spin is most conveniently accounted for by introduc- 
ing the vector c = ((c), (s))"^ as well as a current vector 
/ = {I,Is)^ carrying both the charge and spin on the 
dot and the charge and spin currents into the leads as 
components. The spin on the dot and the spin current 
are defined as the difference of the occupation or current 
for the two spin channels, respectively, (s) = (t) — (i) 
and /, = L 



Pumping current. With the abbreviations introduced 
above, we ca n g ive the expressions for the adiabatic 
charge currenlpS and the adiabatic spin currenlP^l in a 
compact form 



/>.°)(i) = eA(5..)|#°^ 



(46) 



where the components of the instantaneous c are ob- 
tained through 



/t 



= (-l,0,0,l)P, 
= (0,1,-1,0)P, 



(»,0) 



(47) 
(48) 



and Pj*-*'*^-* is the stationary solution of the generalized 



Master equation, Eq. (11). Importantly, in the presence 



of a magnetic field also a spin current is induced by a 
finite bias or a modulation of gates; this spin current 
in turn will in the following be shown to have an im- 
pact on the charge current noise. The current expression, 
Eq. ( 46 ) , shows that time-dependent changes in spin and 



charge on the dot result in a current flow. The amount of 
the latter is found from A, which is a two-by-two matrix 
containing information about the relaxation behavior of 
charge and spin in the dot in response to the external 
driving. This matrix representation of the prefactor is 
owed to the fact that the charge and spin dynamics are 
coupled quantities, i.e., a change of the spin expectation 
value can induce a charge current and vice versa. The 
structure of this matrix is 



A (gaa) 



1^ iaaa) C (ffcHj) 
C,[-9aa) K{-gaa) 



(49) 
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where the functions k, and C are 

f^{9ao) ^ ^^a{Ta+Tgaiy){ga -I) (50) 

a, (7 

C(5a<T) = ^ X!"^^" (-^ + 5'a<T)55^ ■ (51) 



The functions k, C £^nd Q depend on Ta and ^qo- with 
the definition 



(52) 



where the Fermi function for lead a is given as fa {E) — 
1/ (1 + exp[/?(i; - aeV/2)]). We define Q = F (1 - g^g^) 
which is the product of two new relaxation rates occur- 
ring due to the coupling of charge and spin (see, e.g. 
Ref. [83]), divided by F. Also we define a = L if a = R 
and vice versa. When a takes the role of a variable, we 
have a ^ +1 [a = —1), relating to the respective sub- 
script a = L (a = R) and equally for the spin, where 
a = +1 [a ^ -1) relates to f (i)- 

For the specific case where bias and gate are chosen as 
pumping parameters the result for the time-averaged cur- 
rent is shown in Figs, [s] and |9j in each figure (a) and (c). 
We will discuss the shown results in detail in Sec. IIVBI 
Some features can be observed in the pumped cur- 
rent already in the general case. The off-diagonal en- 
tries C represent the coupling of charge and spin dy- 
namics; therefore these coefficients have different sym- 
metries with respect to the Zeeman splitting A, i.e., 
C(-A) = -C(A) whereas k(-A) = k(A). If there is ei- 
ther zero magnetic field or zero interaction, the dynamics 
of charge and spin decouple. It is then straightforward 
to show that the contributions Q vanish and that the ma- 
trix A contains the relaxation rates of charge and spin 
to the different leads a, which, in first order in F are 
given by A^.^^^^ F (1 ± go) (with g^i; = g^i = 5a) with 



A 



(1) 



, — V„ A^'V . For zero magnetic field, 

c/s L^a c/s.cx ^ ' 

arbitrary Coulomb interaction U we find 



A(A = 0) = - 



0, and 



(53) 







This shows that indeed A carries information about the 
relaxation rates. When the dynamics of charge and spin 
decouple, the pumping current is directly related to the 
charge relaxation rates, while the spin current is directly 
related to the spin relaxation rates. Note however that 
for A = the spin expectation value is zero at all times 
and so is the spin current. 

Moreover, if we consider the case of non-interacting 
electrons, we find, independently of the magnitude of A, 



A(C/ 0) = - 



1 Fl - Fr, / 1 



F 



1 



(54) 



The matrix is hence diagonal, and the charge and spin 
dynamics can be considered uncoupled, even with a fi- 
nite magnetic field. This is because in the absence of 
Coulomb interaction the particles with spin f and spin 
I form two independent particle sectors. The only oc- 
curring relaxation rate for [/ — > is then the tunneling 



rate A 



(1) 

c/s 



Fq. Equivalently to the two independent 



particle sectors one can consider any linear superposition 
of them, as in this case the charge (c) + 1 = (t) + (|) 
and the spin (s) = (f) — (J,). An important consequence 
of the vanishing interaction is that the prefactor matrix 
A does not depend on gate and bias anymore. Therefore 
a time-averaged pumped current can be induced only if 
at least one of the F^ depends on time, while the time- 
averaged pumped current vanishes when pumping with 
gate and bias. We will come back to this in Sec. |IVB[ 

The cancellation of the off-diagonal terms in the pref- 
actor matrix for vanishing magnetic field or vanishing 
Coulomb interaction can be observed in the same man- 
ner also in the zero-frequency noise, which we calculate 
in the following. 

Zero-frequency noise. We are interested in the proper- 
ties of the noise indicating the origin of processes leading 
to the appearance or the suppression of a finite time- 
averaged pumped charge. We therefore start by pre- 
senting the results for the instantaneous contributions 
to the zero-frequency noise in the presence of a finite 
bias and an external magnetic field. As already stated 
in Sec. |IID[ the instantaneous noise is found to be the 
time integral of the expression for the stationary zero- 
frequency noise, where all parameters are replaced by 
time-dependent parameters frozen at the integration time 
t. The stationary zero-frequency noise in the sequential 
tunneling regime has been calculated before in Ref. ISTl 
however, analytic expressions were only provided for spe- 
cific limits. Here, we find an analytic expression for the 
instantaneous contribution to the zero-frequency noise in 
a time-dependently driven system 



dt 



A' {gaa 



{go 



(55) 



Importantly, in contrast to the zero-bias case discussed 
before, see Eqs. ( 28 ) and ( 31 1, this expression is no longer 
merely given by the time-averaged conductance. Now it 
contains time-averaged shot noise terms as well; this is 
clear from the first term of Eq. (55 1, which contains the 



instantaneous charge and spin current, /J*'^"*, being non- 
zero only for a finite bias. This term contains the matrix 



(56) 



which has a similar structure as the one occurring in the 
pumping current. The appearing rates are therefore the 
same as in the prefactor of the pumping current, k {gaa) 
and C{gaa), see Eq. (46). Here we observe that these 



rates reappear also in the first term of the instantaneous 
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current-current correlation. The occurrence of rates re- 
lated to charge and spin dynamics in the current noise 
has already been reported for similar systems, see e.g. 



netic field, Eq. (28), has been shown to directly depend 



The second term in Eq. (55) contains the charge and 
spin variance Ac — ((Ac), (As))-^ and for zero bias, it 
reduces to the thermal noise. The charge and spin vari- 
ance enter together with the prefactor vector a that also 
depends on Fq and ga/j, which can be given as 



V {-gaa 



(57) 



For an analytic expression of the function which can 
be expressed in terms related to the relaxation rates of 
the system (like k and C), see Eq. (El ) in the Appendix. 
We have encountered this second term, proportional to 
the local charge and spin fluctuations, Ac, already in a 
simpler form in the previous case. Sec. |III A[ there the 
instantaneous noise in the absence of a bias and a mag- 



on the charge variance and the charge relaxation rate, 
while in the presence of a magnetic field mixing with the 
spin variance occurs. For zero bias, ij*'^'* = 0, only the 
latter term survives (thermal noise). Note however that, 
in the large- bias limit, where shot noise effects are dom- 
inant, both terms contribute equally. Therefore, there 
is no clear separation in thermal and shot noise terms 
in this high-temperature limit with time-dependent driv- 
ing, when representing the instantaneous noise according 
to Eq. ( |55[ ). A similar form for the current noise in the 
time-independent regime has been found in the special 
case for non-interacting electrons, in the zero magnetic 
field and large-bias limit. ^■^ 

Our main interest is here on the first-order Q contribu- 
tion to the noise for arbitrary interaction U and Zeeman 
splitting A. Despite its rather cumbersome analytic form 
it can nonetheless be given in a structured manner, using 
the above introduced matrix notation 



S 



(a,0) 
LR 



dt 



Bl (gacr) It 



dt 



B2 {gaa) li 



7(1,1) . M,0}T 



B3 (5..) J^h 



(58) 



This equation represents a closed form for the pumping 
noise for interacting electrons and an arbitrary choice of 
the driving parameters. We cast the adiabatic noise in 
a form which reflects the structure of the instantaneous 
noise, with the important difference that the adiabatic 
correction to the current and the time derivatives of the 
former contributions appear here. Prefactors of the first 
four terms are collected in the matrices B,;, given in the 



Appendix in Eqs. (|E4|), (|E5|), ^E6^, and ([ETf . These ma- 
trices depend on F^ and ga^, and in particular, B4 con- 
tains the terms proportional to the time derivatives, 
and (jaa- Like A', the B.^ have again the property that 
if charge and spin dynamics are independent, namely for 
U = or B = 0, they become diagonal. 

In addition, we find two terms which contain charge 
and spin variances and their time derivatives, paired with 
the prefactor vectors bi, see Eq. (E8) in the Appendix. 



Here, 62 contains factors proportional to the time deriva- 
tives Fq, and c/aa- 

In analogy to the instantaneous noise, for a vanishing 
magnetic field, A = 0, when spin and charge are indepen- 
dent, the full pumping noise can be expressed in terms of 
quantities related to the charge alone, namely the charge 

on the dot, {c)[^''^\ and its variance, (Ac)^*'^'', the charge 
r(^i) 



current, /j''^'', the charge relaxation rate, xi]a-, and the 
time derivatives of these charge-related functions. 



B. Adiabatic pumping with e(t) and V{t) 

We will in the following use the general results obtained 
above in order to study the pumping current and noise in 
the case where e{t) and V{t) are the pumping parameters. 
We divide the following discussion into two parts. First, 
we recapitulate the discussion of the average pumping 

current, that has been investigated in Refs. I28|29l 
In the second part, we discuss our new results of the 
pumping noise S'lr, and put it into relation with the 
pumping current and the instantaneous current 



1. Pumping current 



The results for the dc pumping current / ' for a fi- 
nite interaction U = 2bk^T are depicted for different pa- 
rameter sets of coupling asymmetry and magnetic field 
in Figs. [S] and |9j (a) and (c). The figures display the 
case of weak pumping amplitudes, i.e., the modulation 
amplitudes of gate and bias are much smaller than k^T. 

The dashed lines in these figures represent the situation 
where either the |0) — \a) or \a) — |2) transition is in 
resonance with one of the electrochemical potentials of 
either the left or right lead. Due to Coulomb blockade 
the charge of the dot is fixed inside the diamonds formed 
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by the dashed hnes. 

The first observation for the pumped charge is that 
there can be only a non-vanishing signal when two reso- 
nance lines cross (whenever two dashed lines meet). The 
reason for this is that only there, the system is effec- 
tively sensitive to two independent pumping parameters, 
namely to the distance between the level position and 
one of the electrochemical potentials, e{t) ± ^V{t). In 
all other regions the system is sensitive to only one of 
these energy differences (namely on the dashed lines in 
between such two points) or not sensitive to the parame- 
ter variation at all (in the regions away from the dashed 
lines). In these latter cases the averaged pumped charge 
is identical to zero. Moreover, the dc pumping current is 
zero for a non-interacting quantum dot, [/ = 0, because 
in that case the charge and spin relaxation rates become 
constant, and consequently there is no difference in the 
loading and unloading process of the pumping current 
(as already stated in the previous Sec. IV A). 



What is of interest are the features of the adiabatically 
pumped current at the crossing points, which depend on 
the magnetic field and the coupling asymmetries. Fig- 
ure [s] (a) shows the pumped charge in the absence of a 
magnetic field, A = 0, and for symmetric coupling to the 
leads, Fl = Fr. The two peaks have opposite sign, which 
was shown to be related to the difference in the degen- 
eracy of neighbouring charge states: for the parameter 
cycle chosen here, the pumped charge is negative, at the 
points where the degeneracy of the ground state is in- 
creased by increasing the charge and the pumped charge 
is positive at the points where the degeneracy decreases 
by increasing the charge. 

Figure 8] (c) shows that a finite pumped charge is 
observed a so at the line crossings at high bias if an addi- 
tional coupling asymmetry to the left and the right lead 
results in an asymmetry between the loading and unload- 
ing process along the pumping cycle, Fl = 0.7F. 

Finally, Fig. [o] (a) and (c) depict the results for the 
adiabatically pumped current in the presence of a finite 
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FIG. 8: Density plot of the pumped charge (a) and (c), 

and the adiabatic noise Sl*^"' (b) and (d), as a function of 
the mean energy level e and the mean bias voltage V, without 
magnetic field, A = 0. Figures (a) and (b) depict the case for 
symmetric tunnel coupling, Fl ~ Fr = r/2, and both (c) and 
(d) for asymmetric coupling, Fl ~ 0.7T and Fr = 0.3F. The 
underlined dashed grid sketches the dot level resonance lines. 
The other parameters are U = 25kBT, 5e = O.IF, 5V = O.IF, 
Q. = O.IF, ip = -7r/2 and F = CS^bT. 
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FIG. 9: Density plot of the pumped charge (a) and (c), 

and the adiabatic noise 5^'^'''' (b) and (d), as a function of 
the mean energy level e and the mean bias voltage V , with 
finite magnetic field, A = 7.5fcBT. Figures (a) and (b) depict 
the case for symmetric tunnel coupling, Fl = Fr = F/2, and 
both (c) and (d) for asymmetric coupling, Fl = 0.7F and 
Fr = 0.3F. The underlined dashed grid sketches the dot level 
resonance lines. The other parameters are as in Fig. [8] 
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magnetic field, A = 7 .hk-oT, for symmetric and asym- 
metric coupling to the leads. First of all, new features 
are observed stemming from the spin-splitting of the en- 
ergy level. Of these additional features, the signals at 
points 3 and 4, as indicated in Fig. [9](d), are visible only 
when the tunnel coupling is asymmetric, Fl 7^ Fr. 

Strikingly, the signal that was observed at point (1) in 
the absence of a magnetic field, vanishes here. This is due 
to the fact that the finite magnetic field lifts the degen- 
eracy of the singly occupied state, resulting in a constant 
relaxation rate and therefore in a symmetric loading and 
unloading in the low-bias regime. As a consequence, the 
pumping current can be written as a full time derivative, 
and its time average vanishes. To show this effect, we 
approximate Eq. ( 46 1 in the vicinity of point 1 for the 



case that point 1 is well separated from all other contri- 
butions (namely under the assumption that A and U are 
sufficiently larger than U-qT) and find 



e Ft, - Fr d 



(59) 



2 F dt' • 
Since the Fq are constant in time, the dc component 
of /("^'O) is zero. Note however, that the prefactor di- 
rectly shows that for symmetric coupling the symmet- 
ric contribution to the time-resolved pumping current, 
^a.o) ^ i^^a.o) _ ^^a,o)^ at point 1 is zcro for all 

times, whereas for asymmetric coupling it cancels only in 
the time average. This vanishing of the pumped charge 
due to the lifting of the degeneracy of th e ground state 
is important as a spectroscopy tool.'^SEII However, from 
the charge current alone it is impossible to distinguish 
whether the pumping current vanishes on average or for 
all times t in the low bias regime. Importantly, as we will 
show later, the pumping noise enables us to differentiate 
between exactly these two cases. 



symmetric point, (e, V") = {—U 12, 0). These relations can 
be proven at the level of Eqs. (46) and (581.^^1 The same 



symmetry relations as in Eqs. (60) and (61) hold also for 



1) 



instantaneous contributions of current, ' and noise, 

Dlr , respectively. 

The complex structure of the pumping noise can be 
more easily accessed when concentrating on working 
points close to the zero-bias region. We start by consid- 
ering the most simple case of symmetric tunnel coupling 
Fl = Fr and zero magnetic field, A = 0, depicted in 
Fig. [s] (b) . In the low bias regime the pumping noise can 
be approximated as 



(a,0) 
LR 



dt 



(i,l)r(a,0) 



(62) 



provided that U 3> k^T. This means that in this limiting 
case the pumping noise can be fully related to the prod- 
uct of the time-resolved currents which are of relevance 
here, namely the pumping current and the instantaneous 
current. This expression containing current expectation 
values, which are again local in time, reflects the struc- 
ture of the current-current correlations in a very sim- 
ple product form. The noise thus contains information 
about the relative behaviour of these two time-resolved 
currents. 

Remarkably, also for asymmetric tunnel-coupling Fl 7^ 
Fr, this simple formula for the low-bias regime still yields 
a good qualitative approximation to the pumping noise. 
Therefore, we will now make use of Eq. (62) to under- 
stand the shape of the noise signal in detail. For this 
purpose we furthermore write the pumping current in 
terms of a response function. A, as introduced in Ref. UHl 



(63) 



2. Pumping noise 

The result for the pumping noise is shown in Figs. [8] 
and [9] (b) and (d), with the same set of parameters used 
for the discussion of the pumping current. In general we 
see that also the pumping noise occurs in the vicinity of 
the dashed line crossings only, since as explained above, 
these are the only points where the system is sensitive 
to two independent pumping parameters. In contrast to 
the pumped charge, the features of the pumping noise 
display a sign change, whenever the pumped charge is 
finite. 

Furthermore, the pumped charge and pumping noise 
generally have the following different symmetry proper- 
ties with respect to the average gate and bias voltage 



,(a,0) 



-e,-V) = -7(-°^(?,F) 

(a,0)._ 



(60) 
(61) 



that is, /'•°'°'' is antisymmetric and ■S'^r^'' is symmetric 
with respect to the point reflection at the electron-hole 



where we used a vector notation for the two pump- 
ing parameters = {e{t),eV{t)) and the gradi- 
ent with respect to said parameters V"^ = (^d^,^dv)- 
The response function for A = is given by A = 

-(a2-A«)/2AW. 

Now, we explicitly insert the pumping current contri- 
bution into Eq. ( 62 ) , where we expand in small pumping 
amplitudes, 



oia-.Q) _ 8-(i,l)-(a,0) 



ASA- 



(64) 

The area of one pumping cycle in parameter space is 
tSA — e /J" dteV. We insert a third component equal 

to zero in V for a well-defined cross product, where the 
third component of the latter appears here. 

The first term in this equation is simply a product 
of time-averaged currents. Since the instantaneous cur- 
rent always changes sign when reversing the bias and the 
pumping current does not change its sign within any res- 
onance, we can generally state that this first contribution 
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has opposite signs for opposite bias with a node at zero 
bias. 

The different shape of the noise signals for different 
tunnel couplings to the leads, see Fig. |8] (b) and (d), 
thus originates from the behaviour of the second contri- 
bution. Since the gradient of the current, V/^*'^\ points 
mainly along the bias axis and the gradient of the charge, 
y ^g^(«,o) ^ jg mainly directed along changes in the gate, the 
cross product has a constant sign for all values of e and 
V. Hence, for a fixed pumping cycle, only the averaged 
response function A is strongly sensitive to the coupling 
asymmetry. It therefore governs the symmetry behavior 
of the pumping noise at low bias. 

For symmetric coupling the response function A has 
a node at zero bias, making the second contribution of 
Eq. (64) of equal shape as the first. This explains the an- 



tisymmetric behavior (with respect to V) of the pumping 
noise at low bias for symmetric tunnel coupling, shown in 
Fig.[8](b). However, when the tunnel coupling is strongly 
asymmetric, as shown in Fig. |8][c) and (d) for Fl — 0.7T 
and Fr = O.SF, the response function A does not exhibit 
any sign change. Therefore the antisymmetric behavior 
of the noise is lifted, see Fig. [s] (d) . 

Concluding the discussion of the A = case, we con- 
sider the pumping noise signal in the high bias regime. 
Also at high bias, pumping noise occurs whenever the 
pumped charge is finite, as can be seen in Fig.[8](d). The 
approximate Eq. (62 1 is no longer valid in this regime, 
and in order to understand the detailed behavior, the full 
Eq. (58) has to be considered. Nonetheless, this pump- 



ing noise exhibits a similar sign change as in the low-bias 
regime, the orientation of which still reflects the respec- 
tive sign of the time-averaged pumping current and the 
instantaneous current (we will describe this effect in more 
detail in the following A case). 

Strikingly, when considering a finite magnetic field, the 
pumping noise can persist even when the pumped charge 
vanishes. In this particular case the pumping noise does 
not exhibit a sign change. This effect can be observed in 
the case of asymmetric coupling in the low bias regime, 
see Figs. |9] (c) and (d). In order to reproduce this re- 
markable feature, we extend our discussion of the pump- 
ing noise in the low-bias regime to the case of a finite 



magnetic field, A 7^ 0. Importantly, Eq. (62) is a valid 



quantitative approximation for finite A independently of 
the tunnel coupling asymmetry, given that ?7, A 3> k^T. 

Since we are here interested in a situation in which 
the average pumped current vanishes, the first term of 
Eq. (64) never contributes and the full behavior can be 



understood from the time-averaged response function. In 
the limit C/, A 3> fceT and for low bias it is given by 



A(A) = 



Fl-Fr 
2F 



(65) 



see also Eq. ( 59 1 . This response function is independent 



of gate and bias, which explains the missing sign change. 
Furthermore, A(A) takes a constant value different from 
zero for a finite coupling asymmetry. For Fl 7^ Fr, the 



time-resolved pumped current is different from zero, even 
though its average vanishes. This explains the finite con- 
tribution to Eq. (|62|). 



However, when the coupling is symmetric, the response 
function A(A) is exactly zero and so is the pumping cur- 
rent at every instant of time. Therefore also the pumping 
noise vanishes, see Fig. [9] (b). This shows that with the 
help of the pumping noise - which, as we want to stress, 
is itself a time averaged quantity - one can distinguish at 
low bias whether the pumping current is zero at all times 
or whether it is only its time average which vanishes. 

Finally, we also remark on the high-bias pumping noise 
occurring in the presence of a finite magnetic field, as 
shown in Figs. [9] (b) and (d). The signals, such as points 
2 to 4, go along with a finite pumped charge and hence 
exhibit the sign change observed and discussed before. 
As mentioned already for the high bias noise signals in 
the case where A = 0, the orientation of the nodes in 
the pumping noise reflect the respective directions of in- 
stantaneous current and pumped charge. Namely, when 
going away from the electron-hole symmetric point, the 
pumping noise cross-correlations change from positive to 
negative (negative to positive) when pump and bias work 
in the opposite (same) direction. This fact also manifests 
in the symmetry relations, Eqs. (60) and (61), derived 
earlier. 



V. CONCLUSION 

In this paper, we developed a formalism for the calcu- 
lation of the zero-frequency pumping noise in the adia- 
batic regime in the presence of strong Coulomb interac- 
tion and a non-equilibrium due to an arbitrary externally 
applied bias voltage. We found analytic expressions for 
the pumping noise containing expectation values of the 
dot occupation and the current, which are local in time. 
This allows us to individuate contributions originating 
from thermal noise, shot noise and pumping noise. 

In a first step we applied the developed formalism to 
the case of pure pumping, in absence of an external bias 
voltage, up to second order in the tunnel coupling F. In 
zeroth order in the driving frequency, the zero-frequency 
noise fulfills a time-averaged version of the fluctuation- 
dissipation theorem. We find that for the correction 
in first order in the driving frequency, the fluctuation- 
dissipation theorem breaks down, uniquely due to the 
nonvanishing Coulomb interaction. This is true already 
in first order in the tunnel coupling F. We study the 
characteristic properties of the pumping noise based on 
an analysis of the adiabatic correction of the Fano factor. 
We find that it exhibits information about the coupling 
asymmetry to the leads and has a distinct feature at the 
electron-hole symmetric point. Interestingly, the adia- 
batic correction of the Fano factor is insensitive to the 
specific choice of pumping parameters. 

In the second part of this manuscript we addressed the 
pumping noise in the presence of a finite - possibly large 
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- bias. In the specific case of pumping with gate and 
bias voltage as time-dependent parameters, we find that 
there can be pumping noise in the absence of pumped 
charge. The appearance of this additional noise signal 
can be used to identify whether the time-resolved pump- 
ing current vanishes or whether it averages out after one 
period. More generally the pumping noise reveals the 
respective direction of pumping current and the current 
induced by a bias. 



and the adiabatic correction (first order in fJ) can be 
evaluated as 

-t-A^"^) (i, z) {t, z) + d,A^'^ {t, z) B(') {t, z) . 

Eventually, for all expressions we encounter, the relevant 
limit is z — > 0+. Nonetheless, we keep z finite at this 
stage, because we will in some cases need to deal with 
derivatives with respect to the Laplace variable z. 
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Appendix A: Adiabatic expansion in Laplace space 

In this section we want to introduce a general scheme 
for an adiabatic expansion, as derived in Ref. ,80,. It 
turns out to be suitable to change to Laplace space. In a 
general case of two functions that depend on two times, 
A (t, t') and B {t, t') we can write the Laplace transform 
of their convolution {AoB){t,f) = J^, dtiA(t,ti)B(ti,f) 
as 

{AoB){t,z)^ f dt'e'<*-'"l f dtiA{tM)B{h,t') 
> Jt' 

dtie-'-^^-^^^ A{t,ti) B iti, z) 

) 

^e"''^''^^ A{t,z)B[t,z) (Al) 

where, as before, the Laplace transform is defined as 
A{t,z) = /^^di'e^(*'-*U(t,i') and analogously for B. 
The operator df is the derivative with respect to z act- 
ing on A and the time derivative of B, respectively. 



From the first to second line in Eq. (|Alj) we swapped 
the integral limits, and from the second to third line, we 
performed a Taylor expansion of B around time t, i.e., 

B{h,z)^Y.n=,hAt^-i^T9tB{t,z). 

The adiabatic approximation takes into account only 
terms up to first order in fi, ergo only first-order 
derivatives of dt are considered. Thus, as long as 
^^^(t, z)dtB{t, z) is sufficiently smalPSl 

[A o B) (t, z) « (1 + d^df) A (t, z) B (i, z) . (A2) 

If each of the objects has an adiabatic expansion of its 
own A{t,z) w A(') {t,z) + A^°-'> (t,z), a consistent ex- 
pansion of the convolution can be given as follows. The 
instantaneous part (zeroth order in Vt) is 



{{A o B) {t, = A« (i, z) B« (t, z) 



(A3) 



In Ref. [Sni this way of writing a general convolution in 
time space was also used for a handy formulation of the 
adiabatic expansion of the kinetic equation. In this case 
the starting point is the equation 

dt'W(t,t')P{t') (A5) 



dP{t) 



dt 



—00 
t 



n=0 



dt'wit,t')Y,^it' 



tY—p{t) 



This expansion is readily identified as the formal series 



dPit) 
dt 



.97 di 



Wt{z)P{t) 



(A6) 



and can serve as the starting point for a frequency expan- 
sion of the kinetic equation in arbitrary order. This equa- 



tion is similar to Eq. (Al) except for the fact that P{t') 
is a function of a single time and, therefore, its Laplace 
transform does not occur in the final result, Eq. ( |A6[ ). 
As pointed out in Ref. EOl this can however be treated in 
complete analogy, when artificially introducing the two- 
time function P{ti,t') = P{ti)6{ti — t' — 0+) and writing 
the kinetic equation as 

dP{t) 



{WoP){t,z) 



dt 

dt'e-'i'^*') I dtiW{t,ti)P{ti)S{ti 
-00 J t' 

t 

dt'W{t,t')P{t') 

— 00 



(A7) 

0+) 

(A8) 



Here the Laplace transform of P{ti, t') is easily found to 
be 



P{t,z) 



dt' e- P{t)5{t ^ t' -0+) 



= Pit) 



(A9) 



In the following we treat functions of a single time when 
occurring in the previously introduced brackets in this 
way, such that 



{(M^oP)(i,z)}W ^ wriz)Pt 



{(iyoP)(t,z)}(") = Wl''{z)P^ 



(0 
t 

M 
t 



(AlO) 
(All) 



+ wr{z)pr+dM''{^)Pt 
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Appendix B: Adiabatic expansion of the 
zero-frequency noise 

Here we want to evaluate the adiabatic correction to 
all terms that contribute to the noise. The instanta- 
neous contribution corresponds to the results obtained 



in Ref. 1741811 where in our case all parameters depend 
parametrically on time. Here we concentrate on the con- 
tributions due to the adiabatic corrections. In order to 



be able to refer to the noise contributions in Eq. ( 10 ) in 
a simple fashion let us name them 



s% it) 



2 T 

e e 



e^e^ 



dti 



dU 



2 T 

-e e 



2 T 

e 



dti / dt 

t pt 
dti / dt2 

to Jti 



dtsWi^ (t, h) H {h,t2) Wi, {t2,h) P its) + (a o /3) 



2 / dhWiJt,h)P{h)(^e'WiAt2,t3)P{t3) + ia^ P) 

to 



dhWi^ (t,ti)P{ti)<E>e^Wi, {t2,h)P{t3) + {a^ (3) 



(Bl) 
(B2) 
(B3) 
(B4) 



r 



where both S^^ and S^^ stem from the term containing 
the product of current expectation values. To obtain 
the two contributions, we simply split the middle time 
integral within the inverval [to,i] into one from tg to ti 
and one from to to t. This separation is used because S^^ 

is of the same structure as Sj^^, and can be treated on 
the same footing. Then we make use of the fact that for 
1^1 ~ ^o| ~^ oo, the propagator takes the form H (ii, to) = 
P (ti) . This replacement we do in both S^p (t) and 

The adiabatic expansion of the first term S^p is done 
in a straightforward manner, by performing the expan- 
sion as in Eq. (A4). Note that in the end, the adiabatic 



expansion is directly analogous to the one of the kinetic 
equation, see Eq. (12 1 



(a) 



2 T 

e e 



n.t^t 



(a) 



<ltPt 

+ (a O /3) 



(0 



Ia.I«,t^^ t 



(B5) 



its adiabatic expansion as 

(a) 



{S^p{t) + S%(t)]\ 



{WjU}f {W,,P} 



(a) 
t 



+ {a^ [3) (B6) 

With the help of the rules for the evaluation of the brack- 
ets we can successively evaluate all contributing objects. 

The adiabatic expansion of the remaining term, 
S^p (t), can be treated in an analogous way only after 
some transformations. By exchanging the first two of 
the time integrals in Eq. (B4) and setting to — > — oo, we 
write 



CO 
t2 



dt2 ee^ dtiWi^{t,ti)P{ti)' 



ee 



dt3Wi^{t2,t3)P{t3) 



(B7) 



We identify in the second line of this equation the ex- 
pression for the current Ij3{t2)- Furthermore, we define 
the two-time object 



The term can be combined with the term S^p, in- 

e''", see 



troducing the object H (t, t') = H {t, t') - P [t) 
Eq. (18 1. The two terms that are combined here are de- 
picted in Figs. [3] (b) and (c) and contain all those terms 
in which the times at which the two currents are taken do 
not lie within the same non-irreducible part of the con- 
tour. We discuss the properties of the resulting object 
n(t,i') in more detail in Appendix [d| The outcome of 
this combination is a convolution of four objects. Thus, it 
can be handled analogously to the term S\ . We perform 



Ia{t, t2 



dtiWiJt,ti)P{ti) , (B8) 



which is related to the current Ia{t) via Ia{t) = Ia{t,t). 
With this definition we are able to write the term S^p (t) 
as the convolution 



W = - / dhL{tM)h{t2) 

J —00 

= -e^'^^'Ut,z)Ip{t) 



(B9) 
(BIO) 
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The only step necessary to fully evaluate this expression 
is the calculation of the Laplace transform of the object 
L{t,t2). We find 



dt,e-^(*-*=) ee^ 



t2 



dti 



= -z(*-tl) 



-[/,(<, z = 0)-/„(i,z) 



dhWi^{tM)P[h) 
(Bll) 



This allows us to perform an expansion of Ia{t, z) in the 
Laplace frequency z. Of this expansion we need the con- 
tribution in zeroth and in first order in z in order to 
evaluate the expansion of the zero-frequency noise in the 
driving frequency VL. We find 



(B12) 



dJait,z) 



-ee^ ^d'^Wi^^tPt (B13) 



The evaluation of S'f| (t) in zeroth and first order in 
the driving frequency in terms of the before introduced 
brackets 



is now straightforward. 



(0 



(a) 



(B14) 
(B15) 



Appendix C: Diagrammatic rules 

In this section, we summarize the rules to diagram- 
matically evaluate the various contributions to the kernel 
and the current kernels as given in Refs. l27l75l81l and the 
new terms occurring in the pumping noise. The example 



shown in Fig. 10 is a diagram that combines features of 
all rules listed below. The rules for the instantaneous 
kernel W^ ^^ (z) in Laplace space are: 

1. Draw all topologically different diagrams with n 
tunneling lines connecting pairs of vertices. Assign 
reservoir index a, energy uj and spin a to each line. 
Also, assign state index x with the corresponding 
energy E^, (t) to each element of the Keldysh con- 
tour connecting two vertices. Additionally, draw 
an external line from the upper leftmost to the up- 
per rightmost corner of the diagram, carrying an 
external frequency —iz. 

2. For each time segment between two adjacent ver- 
tices (they may lie on the same or on different 
branches of the Keldysh contour) assign the resol- 
vent 1/AE (t) where AE (t) is the difference of all 
backward going minus all forward going energies. 




(1) 




FIG. 10: Upper graph: general scheme of a diagram element, 
where t' < t. The index x stands for the different kinds 
of diagrams, W, Wi^, and Wi^i^. The forward part of the 
Keldysh contour is the upper Une going from left to right, 
the lower line represents the backward part of the contour. 
Lower graph: example for one specific diagram first order in 
r, of type Wi^ii^ . The solid arrow line with index a is repre- 
senting a tunneling line where an electron is tunneling in and 
out of reservoir a. The solid dots represent external vertices 
due to the current operator. The dashed lines represent the 
extra imaginary energy lines that are used for either Laplace 
derivatives or adiabatic corrections of diagrams. 



Each vertex containing a dot operator cf^^ gives 

rise to a matrix element ■* |x) where x (x') 

is the dot state entering (leaving) the vertex with 
respect to the Keldysh contour. For the single-level 
quantum dot in particular, the transitions between 
the doubly occupied state and the singly occupied 
state with an electron with spin down, |2) ^ | J,) 
and I 4-) 1 2), pick up a minus sign, —1, due 
to Pauli's principle. All other transitions do not 
acquire an additional factor. 

Each tunneling line with index a contributes with 
a factor of ^Tafa (w) if the line is going backward 
with respect to the closed time path and a factor 
of j^Fq/q (— w) if it is going forward. 



5. Each diagram takes up a prefactor of 
{—i) (—1)'' (—1)*^, where b is the number of 
electron operators (due to internal, namely tunnel- 
ing, vertices) on the backward Keldysh contour, 
and c is the number of crossings of tunneling lines. 

6. Sum over all diagrams that contribute to the same 
kernel element. 

Next, we want to give the additional rules for the 
blocks containing one or two current operators, Wj^\ (z) 



and wlX^,{z). 



The diagrams contributing to Wi^ 
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and Wi^j^ can be directly derived from the rules for 
the kernel W. This direct relation arises due to 
the fact that the current operator, defined as /q, = 
['ya{t)d,l.Caka — h-c] is of analogous structure as 
the tunneling part of the Hamiltonian, see Eq. ([5]), which 
give rise to the tunneling vertices in W. The replacement 
of a tunnel vertex by a current operator (external vertex) 
results in factors ±1 with respect to rule 5, depending on 
the position of la on the contour. Furthermore, an addi- 
tional sign occurs depending on whether a tunneling line 
is incoming or outgoing, due to the shape of the current 
operator. 

The resulting additional rules are set thus that all pos- 
sible replacements are taken into account, i.e., the current 
operator has no fixed position within the diagrams. This 
is possible also for the current operators that are origi- 
nally fixed at time t (see Figs. [2] and [s]), because we take 
the trace e""" on all objects Wj^, Wj^j^ with a fixed cur- 
rent operator. This trace cancels all contributions where 
the external vertex due to /q, is not the first vertex in the 
diagrams, enabling us to neglect the position of la- The 
rules are: 

7. Assign a factor -|-1 for each external vertex on the 
upper (lower) branch of the Keldysh contour which 
describes tunneling of an electron into (out of) lead 
a, and —1 in the other two cases. 

8. Sum up all the diagrams, taking into account the 
above introduced factors ±1 for each possibility to 
replace one (two for W/^/^) tunneling vertices by 
current vertices. For two external vertices, multiply 
with a factor 1/2. 

Finally we give the rules for the corrections of the di- 
agrams in first order in the driving frequency il. In the 
model considered here, there is no need to differ between 
adiabatic corrections of the kernel W and of the objects 
Wi^ as well as Wi^i^. We do on the other hand sepa- 
rate the two cases of adiabatic corrections due to a time- 
dependent tunnel coupling and due to a time-depedent 
energy level .'^^l 

The additional rules for adiabatic corrections due 
to the time dependence of F are: 

9. Add to all diagrams needed for (above) addi- 
tional external frequency lines between any vertex 
ti and the right corner of the diagram, and assign 
to them the energy —izi. 

10. Follow the rules ([T]) to (|6]) taking into account the 
extra lines. 

11. Perform a first derivative with respect to and 
multiply by a factor of ^Fq {t) /^a (t)- Sum all the 
contributions obtained in this way. 

The following are the additional rules for adiabatic cor- 
rections VF*^"^ due to the time dependence of e. 



12. In addition to the external frequency lines added 
according to rule ([9]), put one more external fre- 
quency line from the left corner of the diagram with 
no vertex to the right corner. 

13. Follow the rules ([l]) to (|6| taking into account the 
extra lines. 

14. Perform a second derivative with respect to Zi and 

multiply by -| (^E^ (t) -E^, (t)), where x (x') is 

the dot state entering (leaving) the vertex of the 
external frequency line at U with respect to the 
Keldysh contour. The term E^ (t) {E^> (t)) is omit- 
ted if the segment associated with E^ (t) [E^i (t)) 
does not belong to the diagram. Sum all the con- 
tributions obtained in this way. 

For the purposes of this publication, currents and cur- 
rent noise were calculated in the sequential tunneling 
regime and in the next order in an expansion in F. In 
the sequential tunneling regime, no adiabatic corrections 
to the kernels contribute. Therefore in the present paper 
no adiabatic correction rules for a time-dependent bias 
voltage V are needed. The rules for such a general case 
can be found, e.g., in Ref. [501 For the next-order cal- 
culation, it is furthermore enough to consider adiabatic 
corrections of diagrams up to first order in F. In this 
case, it is possible to formulate a simplification for the 
adiabatic correction in lowest order in the tunnel cou- 
pling, by relating it to the instantaneous kernel. Hence 
we can write down an additional simplified rule: 

15. For the lowest-order term in F simply compute 
IF("-i) (z) = \dM^^^ {A 

Note that this does not hold in general.!^ 



Appendix D: Derivation of 11 
1. Starting point 

In all noise expressions the propagator appears in an 
expression with a finite support, n(t,t') — n(i,t') — 
P(t) ® ^ ■ That is, for the time difference i — i' — ^ oo it 
converges to zero. In Laplace space we can write it as 

n(i,z) = n(i,z) - iF(i) (gje"^. (Dl) 

In order to find a relation to explicitly calculate 11, we 
make a convolution with the kernel W from the left-hand 
side, resulting in 

(VF o n) (t, z)^{yv o n) (t, z)-{w o p) {t, z) (D2) 



The convolutions in the above Eq. (D2) are evaluated 



according to Eq. ( Al I. Since W has no inverse, we need 



an additional condition to fully determine 11, which is 
e"^n = 0. To evaluate the convolution of the kernel W 
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and the propagator H, we have to make use of the Dyson 
equation, relating propagator and kernel 



Adiabatic contribution 



In this section we sketch the derivation of the adiabatic 



n {t, t') = 1+ / dti / dt2W {ti,t2) n {t2, t') (D3) contribution of H. We start again from Eq. (|D2|), where 



where W takes the role of a self-energy. We then trans- 
form Eq. ( D3 ) to Laplace space and obtain 



n(t,z) 



(loVKoH) (i,z) 



(D4) 



where the last term can again be evaluated according to 



equations we need in t 



Eq. (Al). Equations (D2) and (D4| form the full set of 



le following to derive the relations 



that determine the different contributions of 11. 



2. Instantaneous contributions 

Having all ingredients in Laplace space at hand we can 
derive the instantaneous H quite easily. We take the in- 



stantaneous limit of Eq. (D2), i.e., all time derivatives 



vanish and the exponentials carrying these time deriva- 
tives can be approximated by 1. Also, the kernel W as 
well as the density matrix P are approximated to instan- 
taneous order (i.e., zero order in the driving frequency 
17). By taking in a next step the limit z — > 0^ we obtain 



W^n^ = M^/*^n« - hm - 



{i) 



(D5) 



In the last line we used for one the outcome of the , 



limit of the instantaneous Dyson equation, Eq. ( D4 ) 



lim ^n**) (t, z) - 1 = -\. (D6) 



>o+ 



Also, we used the result of the instantaneous master 
equation, = I^t^'^Pj^'^ _ 

The Laplace derivative of 11 can be derived analogously 
by taking again the instantaneous limit of Eq. (D2|. 



However, now we first differentiate it with respect to z 
and then set z to zero 



lim dz 

2-^0 + 



dwl 



zd^W^ 



(D7) 



r(0 



Next, we use the relation H;"' = dW^^Iif^ + W^'^dl^^ 



which is a consequence of the Dyson equation, Eq. (D4), 
in the instantaneous z — ^ O"*" limit. We find 



.(«)tt(*) 



(D8) 



Thus, Eqs. (D5| and (D8) determine the instantaneous 



contributions of H and 9H. 



we e xpa nd all convolutions adiabatically, according to 
Eq. (A4). We then go to the z — >■ 0+ limit of for which 
we get 



[wn] 



(a) 



1 



= {WIi}f' - lim - {WPYf » e^ 



2^0+ z 



{dWP) 



(a) 



(D9) 



Likewise, we perform an adiabatic expansion of the 
Dyson equation for the z — > 0+ limit (note that 



lim. 



H^'^) (t, z) = P, 



_ p(a) 



P 



(a) 



> e^ = {W^n},^'') - lim -P, 

* 2^0 z ' 



1 . 



)e--Hj'. (DIO) 



Inserting Eq. (DIO I into Eq. (D9), and using the adi- 



abatic expansion of the generalised master equation. 



P 



[i) _ 



{WP} 



{M^H}J"^ = P^^^ + H^ ' - {dWP} ^ e^ (Dll) 

which is the desired relation for the adiabatic correction 
to H. 



(a) 



we arrive at 



(0 



Appendix E: Analytic expressions for the noise 
including finite bias and magnetic field 

In this appendix we give the explicit expressions of 
the prefactor matrices and vectors that appear in the 
noise expressions in Sec. |IV A[ In particular, there is the 
prefactor vector a = [ly'igaa-), ^{—gaa-)] in front of Anj*'"-* 



in the instantaneous noise, see Eq. ( [55^ . The function ly 
is found to be 



1 

Q 



gaho 



TlLr, 



go' - gw 



2rh^ 



(El) 



with the additional definitions 



Tn 



^gaa 
LgatT 



(E2) 
(E3) 



Note that for zero magnetic field, we can drop the index a 
and the auxiliary functions ha [h'^) reduce to the charge 
(spin) relaxation rates with respect to lead a. Next, we 
list the matrices appearing in the analytic expression for 
the pumping noise in Eq. ( 58 ) . They can be written as 



Bi = 7^ 



(rVagaa 

^ gaagaa gw ~q' g cirj gaah 



Q 





gaa 







r gacrga 



(E4) 
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Note that for certain cases, A ^ diverges. But because can be seen in Eq. (46 ), the term BiiJ"' ' is well-behaved 
Bi is multiplied with which itself contains A, as 



The remaining matrices are 



B 



B, 



3 — 



^ \ail + 2g„+g^g^) l~2g^+g^g^ J 



The 
Y~n2 



w f l-2g„+g^gi ct (1 + 2g^ + g^.g;) \ 
2 ^ CT (1 + 2.g<, + g-^g^) I - 2g^ + g^g^ J 



-i + .gjf o-(i + 35^) 

Q3 ya{l+g^) -1+ga 



(E5) 
(E6) 

(E7) 



All that is left now are the two prefactor vectors in the 
pumping noise. The one which is preceeding An[^'^^ can 
be given as 



Li R 
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2Q 



r2 
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2r 
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-hgcr 

'^gaK 
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gLagRa 



2Q2 

rg^ / Q 
2Q2 3Q + 4r 
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(E8) 



The second vector, 62, is too big to be displayed in this 
paper. It is of similar shape as 61, but contains terms that 
are proportional to the time derivatives Fq and cjaa- We 
want to stress however, that all results were obtained an- 
alytically, hence there are no additional approximations 
apart from the adiabatic expansion made in Sec. |II D| 
and the lowest-order expansion in the tunnel coupling in 
Sec. mi 
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For the general case where A is the convolution of several 
objects, dA is to be understood as the frequency derivative 



of the convolution at z = 0. 

We would like to remark that in the limit of (/ = the 
two resonant signals of •S'l'^"' simply combine to one sig- 
nal with twice the amplitude. This behaviour differs from 
the instantaneous noise which - as we have pointed out 
in Sec. [Ill A] - gets effectively enhanced in the interacting 
case, due to the enhanced relaxation rate. 
Note that for very large amplitudes one needs to choose 
the frequency Q small enough to be within the adiabatic 
validity regime, QSe <^ TkBT. 

The inversion in the point (—(7/2,0) results in the trans- 
formations gaa — > gaa, C^''"' ~(^t'^\ Ac^'"''' — > Ac^''"' 

as well as /J''^' — > — jj''^', hence the entire noise expression 
remains the same and the current changes sign. Addition- 
ally, since a trajectory in the parameter space of e{t), V{t) 
keeps its orientation under this transformation, the time 
integral does not change the symmetry. 



